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Abstract 

New solutions to the abelian U(l) Higgs model, corresponding to vortices of integer 
and half-integer winding number bound onto the edges of domain walls and possibly 
surrounded by annular current flows, are described, based on a fine-grained analysis of the 
topology of such configurations in spacetime. The existence of these states, which saturate 
BPS bounds in specific limits and are quite reminiscent of D-branes and membranes 
in general, could have interesting and some important consequences in a wide range 
of physical contexts. For instance, they raise the possibility that for some regimes of 
couplings the usual vortex of unit winding number would split into two vortices each 
of one-half winding number bound by a domain wall. A similar approach may also be 
relevant to other known topological states of field theory. 
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1 Introduction. 



The pivotal role played by what is now called the abelian U(l) Higgs model in the progress 
and the stimulus for new ideas in many areas of physics over the second half of the twentieth 
century need not be emphasized Q. Ranging from phase transitions in condensed matter sys- 
tems displaying quantum coherence phenomena such as superconductivity and superfluidity, 
through elementary particle physics with its issues of spontaneous and dynamical symmetry 
breaking, the origin of mass, colour confinement and the dual Meissner effect, to cosmology 
and the evolution of the universe with the formation of textures of different dimensionalities 
all possibly within some inflationary scenario, and finally even M-theory to the extent that 
magnetic vortex solutions provided the first example of so-called BPS states which are so 
crucial to the duality web of superstring theories, the abelian Higgs model has proved to be 
sort of a precious treasure chest for tools in our quest of the secrets of the entire physical 
universe, cast away that we are on our little "speckle of dust". With Dirac's monopole, mag- 
netic vortices were also the first known states bearing witness to the wide-ranging role played 
by the topology of field configurations in classical and quantum field theories alike, thereby 
opening up yet another whole beautiful toolbox of mathematical physics. 

In this latter respect, it seems to be de facto a well accepted working assumption — albeit 
not established in any strict mathematical sense — that all topologically non trivial solutions 
to the Landau-Ginzburg-Higgs (LGH) equations are magnetic vortices of integer winding 
number in the order parameter represented by a spacetime dependent complex scalar field, 

^(x) =/(x)e^^(^). (1) 

Indeed, in a space of infinite extent, finite energy configurations must necessarily be such 
that approaches its non vanishing constant vacuum expectation value at infinity, while 

single- valuedness of ^(x) itself then allows only for a phase dependency of integer winding 
number L at infinity, e~^^'^ {(j) being the angular direction in a plane locally transverse to 
the vortex). Continuity of ^(x) throughout space then implies (when L ^ 0) that iplx) itself 
must vanish along at least one curve of dimension one at a finite distance which determines 
the location in three dimensional space of at least one vortex. Since such configurations 
carry a non vanishing magnetic fiux measured at infinity which is directly proportional to 
the winding number L, they indeed correspond to magnetic vortices. In particular, the 
L = ±1 solution is the celebrated Abrikosov (anti)vortex of superconductors described by 
the Landau-Ginzburg equations!^], or equivalently the Nielsen-Olesen vortex of the abelian 
U(l) Higgs model whose low energy effective relativistic dynamics is that of the bosonic 
string 1^]. Higher winding number solutions with a single vortex are known as giant vortex 
states in superconductivity, and play an important role in the magnetization properties of 
mesoscopic superconductors |^, ^, 0- 

The above topological argument thus provides the classifying scheme for all possible 
vortex solutions (of finite energy) to the LGH equations in two (infinite flat) dimensions, in 
terms of the then necessarily integer winding number of the map from the circle at planar 
infinity onto the U(l) gauge group of phase transformations of the scalar field 'ip{x). Indeed, 
under both the assumptions of finite energy and of an arbitrary collection of discrete zeroes 
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of arbitrary positive integer degree in the order parameter ip{x), it has been shownQ that 
for a specific critical value of the scalar field self-coupling Aq > 0, a critical value which for 
definiteness we shall refer to as Aq = Ac = Ai (the second notation being motivated by the 
relation to vortices of integer winding number), all solutions to these equations are given 
in terms of a collection of such magnetic vortices each of whose positive or negative integer 
winding number equals, in absolute value, the degree of the zero in ip{x) associated to that 
vortex. Furthermore, still precisely for that critical value Ac of the scalar self-coupling, the 
value for the energy of each of these solutions saturates the Bomogol'nyi-Prasad-Sommerfeld 
(BPS) lower bound|P |l^ proportional to its total integer winding number, showing that 
they in fact satisfy first-order differential equations from which their second-order equations 
of motion follow, a tell-tale sign for some underlying supersymmetry. 

The stability of these magnetic vortices with respect to fluctuations within the same 
topological classes of configurations has also been studied for arbitrary values of the scalar 
self-coupling [^, For a coupling less than the critical one, Aq < Ac, and a specified integer 
winding number L, all such giant vortex states have been shown to be stable, thus suggesting 
that an arbitrary collection of vortices of total winding number L would then collapse into a 
single giant vortex with that winding number L. In contradistinction, for a coupling larger 
than the critical one, Aq > Ac, only the L = ±1 fundamental vortices are stable, while 
giant vortices with \L\ > 2 fall apart into a collection of \L\ individual fundamental vortices 
each of winding number sign(L). Finally specifically at the critical coupling, Aq = Ac, all 
configurations corresponding to all possible same sign integer partitions of the total winding 
number L are of equal energy, since they all saturate the same BPS lower bound irrespective 
of the vortices relative positions 14, 15 1. The physical understanding of these 
properties stems from a subtle interplay between, on the one hand, the repulsive (resp. 
attractive) magnetic force acting between vortices whose winding numbers are of the same 
(resp. opposite) sign — much like the force between parallel or antiparallel magnetic dipoles — 
and, on the other hand, the attractive scalar force related to the tendency of the system to 
relax to configurations whose condensate value |'0(3;)| is as close as possible to the constant non 
vanishing expectation value which minimizes the scalar field potential energy (the ratio of the 
latter form of energy density to the former magnetic energy density is precisely given by the 
scalar self-coupling). In particular exactly at the critical coupling, Aq = Ac, the two types of 
forces balance each other, and vortices do not interact with one another whatever their relative 
positions[15|. Clearly, these very properties also explain why for Type II superconductors, 
whose self-coupling is by definition larger than the critical value, Aq > Ac, vortices tend to 
organize themselves into a triangular lattice of Abrikosov vortices, whose lattice constant is 
a function of the scalar self-coupling, as is beautifully confirmed experimentally! 16, 17|. 

In spite of the elegant physical insight offered by these results, the implicit assumptions 
on which they rely may not be satisfied by all solutions to the LGH equations, including 
those of finite energy. Indeed, it is only in the case of the infinite plane that the restriction to 
configurations of finite energy implies that |V'(3;)| must reach its vacuum expectation value at 
infinity in a continuous manner in all angular directions, and it is this latter fact together with 
single- valuedness of ip{x) which then requires the winding number to be integer. However, as 
soon as planar domains of finite spatial extent are considered, a situation which is of interest 
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not only to actual superconductors but clearly also to other systems with phase transitions or 
topological states confined to finite spatial volumes, the value of on the boundary need 

no longer be the vacuum expectation value, thereby allowing a priori even vanishing values, 
and thus also possibly circumventing the apparent restriction to integer winding numbers 
only. What then becomes of the topological classification above? 

Moreover, another assumption always implicit in the results recalled above|8| is that 
the order parameter ip^x) is such that both the corresponding function /(x) as defined in (|l]) 
is positive everywhere, and that if it does vanish, it does so only at discrete points in the 
plane at which the phase value 6{x) is then not defined. However, even though any complex 
number may always be parametrized as in (|^) in terms of its positive amplitude and its 
phase defined only modulo 2tt, when it comes to a single- valued continuous complex function 
defined on the plane (or any Riemann surface) and parametrized as in (|l|), there appears 
a non trivial correlation between the sign of the real function f{x) — which may possibly 
change — and the lack of single- valuedness inherent to the phase parametrization in terms of 
9{x), this correlation being best highlighted by considering the transport of the single- valued 
order parameter ip{x) along any closed Enite contour. Clearly, this last remark is topological 
by nature, and thus provides the means to address the question raised above with regards to 
a topological classification of solutions in finite planar domains. Note also that this approach 
takes the contour of the circle at planar infinity used in the usual argument, to bring it back 
at any point in the plane at a finite distance, thereby providing a finer-grained tool to assess 
the topological properties of solutions to the LGH equations. 

As a matter of fact, it was recently pointed out||l^ that in finite planar domains, beyond 
the usual vortex configurations, the LGH equations also possess annular vortex solutions of 
finite energy and integer winding number, whose order parameter ip{x) vanishes not only at a 
given point, but also on a series of concentric closed curves surrounding that point as well as 
one another in an almost regular radial pattern and at which the function f{x) does indeed 
alternate in sign|^. Each such successive annulus is related to a non vanishing closed current 
flow always running in the same direction, responsible for part of the total magnetic flux and 
whose kinetic energy contributes to the total energy of the configuration an almost identical 
amount, thereby leading to an infinite value in the case of an annular vortex in the infinite 
plane. Hence, these specific properties of such annular vortices explain why they could not 
be uncovered through the general theorems § restricted to only finite energy solutions in the 
infinite plane and such that the function /(x) always remains positive and vanishes only at 
discrete points with some integer degree. 

Finally, and this very point lies at the hear of the whole matter, the above topological 
argument does not take into account the built-in U(l) local gauge invariance of the LGH 
equations which compounds even further the issue of the lack of single-valuedness of the 
function 9{x) for non zero winding number in correlation with the sign of the function f{x). 
In fact, in the same way that for L ^ 6{x) defines a multicovering of the plane as indicated 
for instance by its asymptotic value 9{x) ~ —L(j), it will turn out that f{x) must be considered 
to define a double sheeted covering of the plane, while nonetheless the order parameter 
ip{x) = /(x)e*^(^'' remains continuous, regular and single-valued throughout spacetime. Thus 

^The dynamical and thermodynamical stability of these annular vortices is an open problem. 



3 



for example, a single giant vortex with L 7^ is associated to a degenerate double covering of 
the plane in which the two sheets meet at a single point — namely, the position of the vortex, 
at which / varies as u'^l, with u measuring the distance to the vortex axis — , while the double 
sheeted covering remains nevertheless regular, continuous and differentiable everywhere, with 
the function / having opposite signs on each sheet. 

This paper provides a topological classification of solutions to the LGH equations, 
which we believe should also prove to be complete since it is based on the fine-grained 
topological consideration mentioned above. Clearly, when accounting for the possibility that 
the real function f(x) may a priori take negative as well as positive values in a continuous 
fashion, beyond the ordinary and annular vortex configurations of integer winding number 
L there also exist vortices of half-integer winding number L such that the order parameter 
function f{x) changes sign an odd number of times when transported around some given 
closed contours surrouding such vortices. More specifically, it will be shown that vortices of 
half-integer winding number, which shall be referred to as half-integer vortices for short, have 
the particularity of always being bound onto the edges of an odd number of two dimensional 
domain walls (when viewed in three dimensions) which end on such vortices, and inside of 
which the order parameter ip{x) also vanishes on a two dimensional surface. The same type 
of configuration may in fact also occur for integer vortices (i.e. whose winding number is 
integer), in which case the number of such domain walls must be even. Usual isolated integer 
vortex configurations are thus a particular case of the latter type, with a vanishing number 
of domains walls ending on each of the vortices. Hence, when viewed in a plane locally 
transverse to any such solution, integer and half-integer vortices are characterized by having 
the order parameter vanish not only at the position of the vortex, but also along an even or 
odd number of continuous lines emanating from the vortex itself, in clear contrast again with 
the assumptions of the usual general theorems. Furthermore, in the case of a finite planar 
domain, such vortices bound onto the edges of domain walls may also be surrounded by an 
annular pattern of successive closed current flows extending up to the boundary, in exactly 
the same manner as already described in the case of isolated integer vortices. Clearly, when 
present, such domain walls are associated to the branch cuts of the double sheeted covering 
of the plane which is defined by the function f{x). 

Being bound onto the edges of domain walls, such vortex configurations also display 
yet another remarkable property. Even though the order parameter ^(x) vanishes inside 
the domain wall, close to the edge onto which the vortex is bound a non vanishing closed 
current flow must quantum tunnel through the domain wall, in order to entirely surround 
the location of the vortex and contribute the required amount of magnetic flux. In the 
context of superconductors, this phenomenon is reminiscent of the Josephson effect for S-I-S 
junctions 1 16, 17]. What is a distinctive feature of the integer and half- integer vortices bound 
onto the edges of domain walls however, is that such a quantum tunnel current appears 
within a same single superconductor simply across a surface of vanishing order parameter 
whose location within the material could be anywhere and may also vary in time. 

Since domain walls are regions of space in which the order parameter |V'(a^)| takes values 
which differ from its vacuum expectation value, they necessarily contribute condensation 
energy to the system. In first approximation and for sufficiently extended domain walls. 
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this contribution is expected to be proportional to their length]^. Indeed, their thickness 
is essentiahy constant since it is governed by a specific parameter of the LGH equations, 
namely the superconducting coherence length ^, or the inverse of the Higgs boson mass Mh 
in a particle physics parlance (see section §). Consequently, in the case of spatially bounded 
domains, such domain walls may extend all the way to the boundary of the domain, and 
lead to finite energy configurations nonetheless. On the other hand in the case of an infinite 
spatial extension of the system, these domain walls must necessarily be of finite extent in 
order to lead to solutions of finite energy. In other words, in such a case, these domain walls 
must themselves end on some integer or half-integer vortex within the infinite volume of the 
sample. 

Such considerations thus lead to the following general picture for all the configurations 
which solve the LGH equations. The basic entities in terms of which these solutions may 
be built up are, on the one hand, 1/2-domain walls, and on the other hand, annular current 
flows. By a 1/2-domain wall, we mean a domain wall which ends on a half-integer vortex of 
winding number L = 1/2 or L = —1/2, namely a 1/2-vortex, but whose other edge may or 
may not be bound to another 1/2-vortex, as the case may be. By an annular current flow, 
we mean a closed current flow carrying no winding number, which may be placed around a 
given vortex but then in combination with other such annular flows in the radially outward 
concentric pattern described previously, and which then extends all the way to the boundary 
of the spatial domain (such configurations are thus excluded in the case of an unbounded 
domain, since their energy would be infinite). Indeed, any possible vortex configuration of the 
types described above may be viewed as being constructed from these basic building blocks, 
in exactly the same way that until now isolated giant vortices with |L| > 2 were used to be 
viewed as a collection of |L| individual Abrikosov-Nielsen-Olesen (ANO) vortices stacked on 
top of one another. Of course, this is also consistent with the point of view of the double 
sheeted covering of the plane advocated in this paper. Such a collection of isolated giant 
vortices with |L| > 2 defines a double covering in which subsets of the points at which the 
two sheets meet have coalesced into single points. In the same spirit, the general solutions 
including half-integer vortices and domain walls may thus be viewed as configurations in 
which some of these points and coalesced subsets of points have been split into branch cuts 
in the double covering of the plane, leading in Ene to the basic entities considered above. 

Thus for example, a L = 3/2 vortex may be obtained from three 1/2-domain walls with 
L = 1/2 which share a common edge, while the 1/2-domain walls themselves could or could 
not be stacked on top of one another. Furthermore, each of these 1/2-domain walls could or 
could not end on another 1/2-vortex with either L = 1/2 or L = —1/2, but they must do so in 
the case of a spatial domain of infinite extent in order to obtain configurations of finite energy. 
Indeed, this is the only constraint which must be imposed on the construction, which implies 
that 1/2-domain walls can never extend up to infinity for unbounded spatial domains, while 
they may extend up to the boundary of finite ones. Since the total winding numbers of such 
bound 1/2-domain walls only take the values L = 0, ±1, the fact that only such basic entities 
are allowed in the case of unbounded spatial domains is in full agreement with the restriction 
to only integer total winding numbers at spatial infinity as implied by the usual topological 

property reminiscent of the close-to-linear quark potential in QCD, as S. Vandoren pointed out. 
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argument recalled earlier. Nevertheless, the assumptions of the usual theorems|g|, 12] are 
evaded since the order parameter then vanishes not only at discrete points in the plane but 
also along segments of finite length, thereby allowing for the existence of these new types of 
configurations built up from bound 1/2-domain walls (annular current flows are excluded for 
unbounded domains). In contradistinction in the case of a planar domain of finite extent, 
the total winding number need not be integer, but could also take half-integer values, since 
a 1/2-domain wall may extend all the way to the boundary of the planar domain without 
being bound to another 1/2- vortex, leading to a finite energy configuration nonetheless. 
Furthermore in this case, such configurations may also be combined with annular current 
flows surrouding vortices and extending up to the boundary of the finite domain. 

The characterization described above of all solutions to the LGH equations, which is a 
direct consequence of the finer-grained topological analysis based on any contour at a finite 
distance in contrast with the usual argument, also raises the issue of the stability and lowest 
energy configuration of bound 1/2-domain walls. Given the fact, recalled previously, that 
magnetic vortices whose winding numbers are of the same (resp. opposite) sign repel (resp. 
attract) each other, while any deviation from its vacuum expectation value in the order 
parameter iplx) induces an attractive scalar force, one should expect that the bound 1/2- 
domain wall binding two 1/2-vortices of opposite winding numbers L = 1/2 and L = —1/2 
would collapse, in its lowest energy state, to a L = configuration without any zero in the 
order parameter. On the other hand, the same issue in the case of bound 1/2-domain walls 
binding two 1/2-vortices with the same winding number L = ±1/2 is far less obviou^. Given 
the natural tension built into such vortex and domain wall states, their energy must increase 
with their curvature or bending, so that the lowest energy configurations are expected to be 
straight vortices and domain walls. Nevertheless, as pointed out above, for sufficiently large 
separation the energy of such states must also grow essentially linearly with the distance 
between the two 1/2-vortices, 1/2-domains behaving very much like rubber bands in such 
a regime because of the condensation energy stored in the domain wall. Hence, it is only 
when the two 1/2-vortices have a non vanishing overlap that the competition between the 
repulsive and attractive forces acting on them could be such that the equilibrium lowest 
energy configuration would be reached for some specific non vanishing separation rather than 
a fully collapsed state. The possible overlap between vortices depends on their size, which is 
governed by yet another parameter of the LGH equations, namely the Meissner penetration 
length A, or the inverse of the gauge boson mass in a particle physics parlance (see 
section |^. In fact, the ratio of this latter length scale to that which determines the width 
of domain walls is directly related to the value of the scalar self-coupling. Consequently, 
one would expect that there might exist some critical value A1/2 of the scalar self-coupling 

^Note that such domain waUs would repel or attract one another as well, according to whether the vortices 
which are bound onto their edges have winding numbers of identical or opposite signs. For instance in the case 
of the 1/ = 3/2 example briefly described previously and assuming that these bound 1/2-domain walls do not 
collapse into Abrikosov vortices, one should expect for this reason that in the infinite domain case the lowest 
energy configuration for the L = 3/2 vortex is obtained when the three 1/2-domain walls are not stacked on 
top of one another, but rather are spread out into a maximally symmetric star-like pattern with an angular 
openin g of 120° between each, the complete state then having a total winding number L = 3 (see however 
section 3.1 ). 
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marking the boundary between, on the one hand, unstable bound 1/2-domain walls which 
would always collapse to L = ±1 ANO vortices whatever the separation between the bound 
1/2- vortices, and, on the other hand, stable bound 1/2-domain walls which would stabilize 
into configurations such that the domain wall is streched to some finite length. In fact, given 
the situation which applies to the stability of integer vortices against fluctuations which are 
also of integer winding number one would expect that it is for values of the scalar self- 
coupling Ao larger than this critical value A1/2 that the latter occurrence would be realized. 
Furthermore for the same reason, it may even be possible^ that the two critical values Ai = Ac 
and Ai/2) thus characterizing the stability properties of isolated integer vortices and of half- 
integer vortices bound by domain walls, should be equal to one another, and thus be equal 
to the value Ac for which the usual isolated integer vortices also saturate the BPS bounds. 

Such considerations thus also raise the possibility that the L = ±1 ANO vortex could 
be unstable against splitting into two L = ±1/2 vortices bound onto the edges of a bound 
1/2-domain wall with L = ±1, and which would stabilize from one another at some specific 
non vanishing distance whose value would depend on that of the scalar self-coupling. Even 
though some arguments and indicative results will be presented concerning this issue, no 
conclusion can be drawn at this stage unfortunately, since it requires a full-fledged dedicated 
numerical analysis which still needs to be completed. Obviously, this specific point, and in 
particular the value for the possible equilibrium distance as a function of the scalar self- 
coupling, is of great interest, for instance, to the physical properties of magnetic vortices in 
superconductors. For example, given a value for that coupling such that the ANO vortex 
would be split, the Abrikosov lattice acquires a fine-grained structure in which the centers 
of all bound 1/2-domain walls still define the usual triangular lattice, but with all individual 
bound 1/2-domain walls then all parallel to one another and aligned with one of the three 
symmetry axes of the fundamental triangular lattice cell, as suggested by symmetry and 
minimal energy considerations. 

The above analysis, based on the fine-grained topological argument involving closed 
contours at any finite distance, has thus unravelled quite an unexpectedly rich classification 
of vortex, domain wall and annular current flow configurations (the latter only in bounded 
spatial domains) of integer and half-integer winding numbers, solving the LGH equations and 
extending the well established isolated vortex solutions. Such states can all be built out of 
some basic entities bound and stacked together, namely the 1/2-domain wall and the annular 
current flow in the case of bounded spatial domains, and the bound 1/2-domain walls in 
unbounded spatial domains. Even though the bound 1/2-domain walls of winding numbers 
L = ±1 may not be stable against collapse into the L = ±1 ANO vortices for some regimes 
of scalar self-coupling, nevertheless the mere existence of these fluctuation modes of bound 
1/2-domain walls with winding numbers L = 0, ±1 solving the LGH equations, must have 
important consequences with regards to the dynamical properties of systems described by the 
LGH equations when coupled to external sources of electromagnetic disturbances, the first 
example in point being again obviously superconductors. 

^However, one should not overlook the fact that the stability analysisjl^ of isolated integer vortices con- 
sidered only fluctuations within the same topology class, thereby not necessarily accounting for the possibility 
of half-integer vortices nor domain walls. 
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The integer and half-integer winding number configurations bound by domain walls 
considered in this paper are also characterized by specific cut structures in the functions f{x) 
and 9{x), and a double sheeted covering of the plane or a finite domain of it (see section 
p.4| ). Such cuts however, have no bearing whatsoever on the physical consistency of these 
configurations, since not only does the order parameter ip{x) remain single- valued, regular 
and continuous throughout spacetime, but also it is only = f'^{x) which is a gauge 

invariant physical observable, measuring in the context of superconductivity for example, 
the local (relative) density of Cooper pairs, and in the particle physics context, the local 
condensation of the scalar field. As will become clearer in the course of the analysis, it is the 
U(l) gauge invariance properties of the system which provide the essential guarantee for the 
physical consistency and existence of half-integer winding number and domain wall solutions 
to the LGH equations. 

In recent years, similar half-integer and fractional winding number topological confi- 
gurations have in fact also been discussed within the framework of models other than the 
abelian U(l) Higgs model, which could possibly have applications in high temperature su- 
perconductors and superfluids|l^, in the QCD confinement problempO|], and for so-called 



Alice stringsj21| in models beyond the Standard Model of particle physics. What distin- 
guishes all these other instances of fractional winding number from the half-integer vortices 
discussed in this paper, is that these other models possess some spontaneously broken non 
abelian internal symmetry, whether global or locally gauged, such that when taken around 
a closed contour the matter fields of these systems end up pointing in a different direction 
in the representation space of their internal symmetry. What is remarkable with the present 
half-integer vortices is that a similar mechanism is at work in spite of having only an abelian 
U(l) gauged symmetry, by carefully keeping track, when taken around closed contours, of 
the possible change of sign in /(x) as well as its lack of single-valuedness in correlation with 
that of the phase parameter 9{x), while the order parameter tp(x) itself remains nonethe- 
less single-valued throughout space. As will become clear hereafter, this particular feature 
is achieved through the U(l) gauge invariance properties of the considered system and the 
correlation which exists between its gauge dependent variables. 

The paper is organized as follows. In the next section, the general equations of the 
abelian U(l) Higgs model will be considered, to lead to the topological classification of in- 
teger and half-integer vortices. Section |3| then particularizes to the two dimensional case, 
namely the plane transverse to such solutions, to discuss specific limits and properties re- 
lated to the possible BPS bounds these solutions could saturate. Single half-vortex solutions 
in finite disks and annuli are then considered in section ^ including the results of some pre- 
liminary numerical simulations. Similar issues are then addressed in section ^ in the case of 
the bound 1/2-domain wall, including some preliminary indicative numerical results as to the 
possible instability of the L = I vortex. Finally, open issues and further outlook onto possi- 
ble consequences and applications of these topological configurations of half-integer winding 
number of the abelian U(l) Higgs model are discussed in the Conclusions. 



2 The Topological Analysis 
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2.1 The LGH equations 



In particle physics units, such that h = 1, c = 1, eo = 1 and /io = l/(eoc^) = 1 — cq and /xq 
being the vacuum electric and magnetic permitivities, respectively — the abelian U(l) Higgs 
model is defined by the following action]^, 

loo) {"I^'"'^'^- + I(Vm + iQA^)i^ |2 -^Ao (l V I' -a')'} , (2) 

with Ffj^i, = Vfj,Ai, — Vu^fi and = d/dx^ (/i = 0, 1,2,3; the notation is reserved for 

another purpose hereafter), while our choice of Minkowski metric signature is (H ). 

Here, q stands for the U(l) charge or gauge coupling of the complex scalar field il^{x), and 
a > for its expectation value with the units of a mass scale which is then also the dimension 
of the field il){x) itself. 

By construction, this system possesses a U(l) local gauge invariance under, 

^'{x) = e'^^'''^^{x) , A'^{x) = A,{x)-^V,x{x), (3) 

x{x) being an arbitrary function which is regular throughout spacetime. This symmetry is 
spontaneously broken through the Higgs potential of dimensionless scalar field self-coupling 
Ao > 0, thereby leading to a massive gauge boson 7 and a massive Higgs scalar h both of 
zero U(l) charge such that 

My = \qa\ , Mh = V2}^a. (4) 



As is already implicit in most of the discussion presented in the Introduction, we shall 
rather work with another choice of units, directly appropriate to an application to supercon- 
ductors. Nevertheless, a translation from one choice of units to the other is straightforward 
enough, for which the rules are given hereafter. Normalizing the scalar field tp — representing 
the Cooper pair wave function — to its vacuum expectation value — namely the square-root of 
the BCS condensate in the bulk in the absence of any magnetic field — , the corresponding 



action is given by|ll8|| 



eoc 



dt 



(fx 



1 



1 / 



{00) 



2 V27rA 



2 r 



(V 



2^ 



Ml* 



-1 



, (5) 



*Note that the specific case of a four dimensional spacetime is considered, but that the results of this 
paper may easily be extended to higher dimensions and thereby lead to domain wall solutions of integer and 
half-integer winding numbers of arbitrary spacelike dimensionalities and bound to the edge of domain walls of 
one more spacelike dimension. Such configurations are very much reminiscent of the D-branes^^ in string and 
M-theory (see also section B.3), thus raising many further issues of interest which are beyond the scope of this 
work. In particular, in the same way that the low energy effective dynamics of the ANO vortex is described 
by the Nambu-Goto actiouQ, the same dynamics for such domain walls and those considered in this paper is 
necessarily also provided by the Nambu-Goto action extended to world- volumes of the appropriate dimension. 
In the same vein, one could also consider possible generalizations to some complex p-form coupled to a real 
g-form with a spontaneously broken local symmetry. 
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where now all quantities are expressed in S.I. units, with of course = {x 
A^^ = {^/c,A) — $ being the electromagnetic scalar potential — and F^y 
the above expression, $o stands for the unit of quantum of flux 



and x' 







2-Kh 



>0, 



ct, 
In 



(6) 



q = — 2e < being of course the Cooper pair electric charge, while A and ^ are the tempera- 
ture dependent magnetic (and electric) penetration and coherence lengths, respectively. Note 
how these two length scales indeed measure the relative contributions to the effective action 
above, and thus also to the associated energy, of the magnetic and condensate energy densi- 
ties. Indeed, the scale ^ weighs the contribution of the gauge covariantized gradient of any 
deviation in ■i/' from its vacuum expectation value against its potential energy, while this to- 
tal condensation energy in turn is measured against the electromagnetic energy contribution 
through the scale A. Note also that one has. 



c 



d <^ 1 9 „ 


dx^ c cdt 



ijk 



d 



d 
dx^ 



-e^^^F^' , i,j,k = 1,2,3, 

where E and B stand of course for the electric and magnetic fields, respectively, and e^^^ is 
the totally antisymmetric Levi-Civita tensor in three dimensions with e^^^ = 

Comparing now the above two actions, the translation between the particle physics 
units and the "superconductor" ones is as follows. 



a <- 



2ttX 



M ^ M V2 Mh r- 



1 



(8) 



where k = \/^ defines the Landau-Ginzburg (LG) parameter of the superconductor. In 
particular, the critical value Ac = Ai for the scalar self-coupling mentioned in the Introduction 
then corresponds to the value with Kc = 1/^/2, which, in the particle physics 

context, thus implies M/^ = M-y corresponding to Aq = Ac = Ai = 1. 

In terms of the "superconductor units", the LGH equations then read as follows. For 
the LGH equation proper, we have the covariantized Landau-Ginzburg equation 



n 



+ ^ (IV-P - i) V' = 0, 



which is thus coupled to the inhomogeneous Maxwell equations^ 



em /i: 



(9) 



(10) 



where the conserved electromagnetic current density 
0, is given by (-0* being the complex conjugate of tp) 

ih 



such that fiJlioi 



em fi 



2gA2 



r 



(11) 



The homogeneous MaxweU equations are of course a consequence of the definition of the field strength 
F^v in terms of the gauge potential y4p, being the corresponding Bianchi identities. 



10 



As a matter of fact, still another choice of normalized units is far more convenient, using 
the canonical length scale provided by A and the canonical magnetic field value provided by 
<l>o/(27rA^). Let us thus introduce the following normalized spacetime coordinates 



ct ^ X 



such that 



V, 



d 



d 



-df, , d^ = {do,d) , do = dr = -^ , d = — 



A 



(12) 



(13) 



while in the electromagnetic sector, the normalized quantities e, b, ffj,u, = (</?, a) and 
J'^ = {J^, J) are defined as follows. 



B = — %b 



S _ $0 
7 ~ 27rA2 



27rA2 



\0 







1 ^0 ^0 



2ttX 

1 ^0 



if , A 



27rA 



J. 



" A 27rA2 ' ^'""''^ ~ A 27rA2 
In terms of these normalized variables, the LGH equation then reads 



(14) 
(15) 



{d^ - ia^) {d'' - ia^') = (l- | ^ V or [(9^ - iipf - {d + iaf^ V = (l - IV'P) 4', 

(16) 

while the inhomogeneous Maxwell equations are 



d'^fuij. = Ju or d ■ e = , d xb — drC = J, 



(17) 



with the components of the Lorentz covariant electromagnetic current density = (J^, J ] 
given by 

(18) 



Ju 



or 



J" 



1 



and satisfying the conservation equation, 

a„J/^ = or drJ^ + d- J = 0, 



(19) 



(20) 



which as always follows from the inhomogeneous Maxwell equations. Moreover, we then also 
have 

ftiv = d^au-dua^ or e = -d^p - drd , b = d x a, (21) 

relations which imply the homogeneous Maxwell equations as the corresponding Bianchi 
identities, 

e^^P^d^fpa = OT dxe + drb = , d-b = 0, (22) 
gMi'po" being the four dimensional Levi-Civita totally antisymmetric tensor with e0i23 = 
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Let us now also introduce the Lorentz covariant current j'^ = whose components 

are defined by dividing the current density {—J^) opposite to the electromagnetic one by the 
square of the spacetime local condensate value IV'I, 

J. = -^ = ^^5Mln(|^)+aM, (23) 

or 



In particular, one then finds 1 23 1 



ftiu = d^ju - dujf^ or e = -drj-df , b = dxj. (25) 

Note that the last of these relations is nothing but the celebrated second London equation, 
which explains the Meissner effect in superconductors. On the other hand, the second of these 
relations is the appropriate Lorentz covariant extension of the celebrated first London equa- 
tion. Indeed, the abelian U(l) Higgs model provides an effective theory of superconductivity 
which is also Lorentz covariant, and in which electric and magnetic fields play roles dual to 
one another under Lorentz boosts. Contrary to the usual LG and London equations which are 
not Lorentz covariant, and do not allow for electric fields in superconductors, the covariant 
formalism must necessarily imply new phenomena in specific regimes where electric fields can 
no longer be ignored. Clearly a covariant formalism implies specific differences in the time 
dependent dynamics of superconductors, but it does so already for static configurations of 
applied electric and magnetic fields, as discussed recently in Ref.|23]. 



Finally, let us give the expression for the free energy of the system in our choice of 
normalized units, when subjected to some external electric and magnetic fields Sext and 6ext- 
One has 



b-h 



'cxt 



2 

+ 



+ \idr- iip)il^ ? + \{d + ia)il) ^ +1^^ (l- | ^ - h'^] ■ 

(26) 

Note that by having subtracted in the integrand the constant term in fi:^/2, we have chosen 
to fix the zero value of the free energy at the normal-superconducting phase transition for 
which ■0 = and e = Cext; b = 6ext- Depending on the consideration we wish to emphasize, 
especially when discussing solutions in the infinite plane, sometimes we shall not include that 
constant term in the evaluation of the total energy of the system. 

The physical interpretation of the overall constant factor multiplying this expression 
for the free energy should be clear, since it corresponds to the magnetic energy of a constant 
magnetic field of value $o/(27rA^) over a volume A^, precisely the two scales with respect to 
which all quantities have been normalized. 

It is also useful to revert to particle physics units given our choice of normalized quanti- 
ties. From the translation rules given previously, it follows that distance scales are measured 
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in units of 1/Mj, magnetic and electric fields in units of M^/\q\, and finally the energy in 
units of 

thus displaying the well-known fact that solitonic solutions to field theory equations have 
mass values which are proportional to the gauge boson mass scale and inversely proportional 
to the squared gauge coupling constant. 



2.2 The polar parametrization 



Given the different expressions above, let us now consider the polar decomposition of the 
scalar field, 

V(x)=/(x)e*''(^), (28) 

where f{x) is thus a real function and 9{x) is an angular variable defined modulo tt rather 
than 27r. Indeed, since the sign of the real function f{x) cannot be constrained to be always 
positive, which, for reasons discussed in the Introduction, would be an unjustifiable restriction 
in the case of a complex function as opposed to a complex number, the arbitrariness in the 
overall sign of f{x) translates into a definition of 9{x) only modulo vr (the reason for this 
arbitrariness will be fully understood in section on basis of the double sheeted covering 
of the plane mentioned in the Introduction). Moreover, the lack of single- valuedness in this 
quantity, which is thus correlated to the choice of an overall sign for the real function f{x), is 
also compounded with the U(l) local gauge freedom inherent to the system, an issue which 
shall be addressed more closely at the end of this section and in the next one. 

The conserved electromagnetic current density J'^, the source term for the inhomoge- 
neous Maxwell equations, is then given by 



J^. = -fj^. or 



2 .-0 



J 



with thus the conservation equation 

{ff) =0 or dr (ff) + d ■ (fj) = 0. 
In particular, the phase variable 9 is determined from the equations (|2^) and (|24 



-Jti + or drt 



89 



a. 



(29) 

(30) 
namely, 
(31) 



Separating the real and imaginary parts of the LGH equation (16), one finds that for / 7^ 
the imaginary part is equivalent to the current conservation equation as befits a U(l) gauge 
invariant theory, while the real part reads, 

,o2^ 



f={r-j'')f-Ai-f)f- 



(32) 



In order to identify an organizational principle within this set of coupled differential 
equations, let us now substitute the covariant London equations (^) into the inhomogeneous 
Maxwell equations. This leads to. 



j = fj + d(drf + d-j 



(33) 
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The expected appearance of the Klein-Gordon operator d^^d^ = [d^ — 9^] in these different 
equations suggests to also establish a similar equation for 9, given the relations in (^T|), leading 
to 



drf + d-j)-{drV + d-a). (34) 



In conclusion, we have thus obtained the following set of equations, which may be 
expressed either in a manifestly covariant form or in terms of each of their components 
separately. The Lorentz covariant time dependent LGH equation is thus 

d>.d>^f = j^ff + K\l-f)f, (35) 



or equivalently the equation given in (32). This LGH equation is coupled to the inhomoge- 
neous Maxwell equation^ 

dud''j, = -fj, + d,{d,f), (36) 



or in component form the equations in (|33D, from which the current conservation equation 
( ^0|) follows. Once this coupled set of equations for / and j'^ = solved (for a given 

choice of boundary conditions to be discussed presently) , the associated electric and magnetic 
fields within the superconductor are determined from the Lorentz covariant London equations 

ffiu = d^ju - dui^ or e = -drj - df , b = dxj, (37) 



from which the homogeneous Maxwell equations (|2^) follow. One then needs to identify 
gauge potentials a'^ = {ip, a) associated to these two fields, such that 

ffj.u = d^au-duafj_ or e = -d(p-dra , b = d x a, (38) 



to also determine the phase variable 9 through the equations given in (|3^) , or the second-order 
equation in (|3^, namely 

d^,9 = -j^ + a^ , a^9^e = -9^j^ + 5X- (39) 

Finally, the free energy of the system may also be expressed solely in terms of the variables 
/, and J, the electric and magnetic fields within the superconductor being obtained from 
the covariant London equations (E^ 



TP - ^ 

^ ~ 2mo 



+ 



+ {drff + [df)' + [f'+?) P + (1 - f)' - i^^l . 

(40) 

The above relations thus determine the system of coupled equations which is to be 
solved, once the relevant boundary conditions specified. Outside the superconductor, the 



^These equations may also be expressed in first-order form as d'^fv^ = — /^jp, provided one uses tlie 
relations between tlie electric and magnetic fields and tlie current j'^ — whidi are defined by tlie 

covariant London equations (|25[). 
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scalar field ■0 vanishes of course, and since the covariant London equations no longer apply, 
the relevant complete set of inhomogeneous and homogeneous Maxwell equations then has to 
be considered on its own, including the possibility of applied fields or sources for these fields 
providing further boundary conditions. At the boundary with the superconductor, matching 
conditions for the electromagnetic field strength f^i, and potential must be enforced. 
Moreover, when considering variations of the action or the free energy of the system with 
respect to the order parameter ^0, one concludes that the current (9 + ia)^l^, namely the 
gauge covariantized space gradient of ip, must have a vanishing normal component at those 
sections of the superconductor boundary which are in contact with an isolating material. 
Separating the real and imaginary parts of that quantity, this condition thus requires that 
both the electromagnetic current J = —f'^j as well as the gradient df must have a vanishing 
component normal to such boundaries. These different conditions would thus appear to 
provide the complete set of necessary boundary conditions, with one exception however, 
to be discussed in the next section. Indeed, the integrated London equations, through the 
associated electromagnetic flux values, provide further constraints of a global character which 
encode the topological vortex structure within the superconductor, which the local form of the 
covariant London equations in (^) cannot account for. These "global boundary conditions" , 
are the last essential conditions which should define unique and gauge invariant solutions to 
the LGH and Maxwell equations above. 

Until now, we have been careful to maintain manifest the Lorentz covariance properties 
of the system throughout. It should be clear that in the limit of time independent configu- 
rations, namely static or stationary ones, the above set of equations reduces to the familiar 
Landau-Ginzburg equations for superconductors solely submitted to magnetic fields (thus 
with e = 0, if = and = 0). However, also for reasons advocated in Ref . ||23[| , when consi- 
dering time dependent configurations or situations involving electric fields as well, the usual 
equations can no longer be applicable, since in particular the usual first London equation, 
drj = —e (rather than the covariant one above, drj + dj^ = —e), implies that superconduc- 
tors cannot sustain electric fields in static configurations, in clear contradiction with Lorentz 
covariance. Another way to realize that the usual LG equations cannot be Lorentz covariant 
is to consider the non relativistic limit in which the parameter c is taken to infinity. Since 
in S.I. units electric fields are measured in units of E/c relative to magnetic fields B, in 
that limit both any time as well as any electric field dependencies decouple from the above 
equations, which then also reduce again to the usual time independent non covariant LG 
equations. It is thus of great interest to explore the physical consequences in relativistic and 
electric regimes of the abelian U(l) Higgs model as the canonical Lorentz covariant extension 
of the usual Landau-Ginzburg effective field theory description of superconductivity |]^. Note 
however that the solutions which are described in this paper are static ones, in the absence of 
any electric field, and thus solve in fact also the usual non covariant LG equations coupled to 
Maxwell's equations for magnetic fields, even though these states in fact define Lorentz co- 
variant configurations solving the LGH equations when subjected to arbitrary Lorentz boosts 
acting from their rest frame. 

The advantage of having separated the system of coupled equations through the orga- 
nization just described is not only that first one needs to solve only for the quantities /, 
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and j, and only then for the phase variable through the knowledge of the electric and mag- 
netic fields and their associated gauge potentials]^. More importantly perhaps, for obvious 
physics reasons the quantities f'^j^ and /^j are in fact gauge invariant variables, directly 
amenable to physical observation, at least in principle, which is a property also shared by 
the electric and magnetic fields e and 6 only. Indeed, U(l) gauge transformations of the field 
variables are given by, 

^' = e'''^ , A'^ = Af,--V^x , a'^ = a^ + d^x, (41) 



so that 



Q 



+ X , f =^ + drX , a =a-dx, (42) 



where xi'T^'^) is thus an arbitrary local function which is regular throughout spacetime. As 
was discussed above, the real function / is defined only up to an overall sign whose choice, 
however, leaves all the above equations invariant as it should, and which indeed does not lead 
to any physical consequence. Nevertheless, this choice of overall sign in / is also correlated 
to a constant shift by zbvr in 9, again without consequences for physical observables. This 
correlated arbitrariness in both / and 6 is then compounded further in the case of 9 because 
of the U(l) gauge symmetries of the system, which induce arbitrary spacetime dependent 
shifts in that variable in correlation with specific transformations in the gauge potential a^. 



In particular, note that the equations ( |3lD or (^) and (3£) display this correlation explicitly 
through their manifest gauge invar iance. 

Consequently, except for the choice of overall sign in the function /, which is without 
physical significance, all the gauge dependency of the system resides in the variables 0, ip 
and a. The remaining variables, namely / — except for its overall sign — , , j, h and e 
are all gauge invariant quantities, whose determination through the above equations is thus 
decoupled from that of the gauge dependent variables. Hence the advantage in using the 
procedure outlined above for the resolution of the coupled system of equations. In fact, as 
will be discussed in the next section, the U(l) local gauge freedom is such that there always 
exists a choice of gauge for which the function 9 is then uniquely determined in terms of 
a specific non regular function which completely encodes the entire topological vortex 
structure of a given configuration, and in terms of which the gauge potential is then 
uniquely constructed from (|39|) as = d^9o + j^. Hence indeed, through an appropriate 
choice of gauge associated to a specific vortex configuration in spacetime, only the LGH and 
inhomogeneous Maxwell equations (32) and (|3^) remain to be solved for the quantities / and 
j/j. — (^j^^j^j subjected to the appropriate boundary conditions. 

To conclude, let us also make the important following point. Even though the currents 
and j are gauge invariant, they are not directly amenable to physical observation, since 
they are so only through the electromagnetic current density {—/"^j^, —f'^j)- Hence, it is only 
at those points in spacetime where / does not vanish that the current is well-defined, 

while at those locations where / does vanish, it may well be that this current possesses 

*As will be discussed in the next section, one should keep in mind however, the subtle connection between 
these two sets of variables which involves the vortex topological structure through the integrated London 
equations. 
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singularities just mild enough to be screened by the vanishing order parameter ip. As we 
shall see, this is precisely what happens not only at the location of vortices, but also on the 
domain walls ending on them. In fact, the expressions in ( p3| ) and (|2^) are already indicative 
of the logarithmic cut structures that may appear in j^, j, f and 6 at the zeroes of the 
order parameter ^ in a way which depends on its winding number at these points. Note 
also that the same remark applies to the singular character of the gauge dependent potential 
at the location of vortices, a well known property typical of topological configurations 
related to gauge fields which is made explicit in the present instance through the appearance 
of the singular function 9q in the determination of a^, while electric and magnetic fields are 
nonetheless well defined throughout spacetime, including at the location of vortices. 

2.3 Configurations of integer and half-integer winding number 

Let us now apply the fine-grained topological analysis described in the Introduction to the 
general framework which has been set-up so far. For this purpose, let us consider any closed 
contour C in spacetime and measure the total electromagnetic flux through any two dimen- 
sional surface S with C as a boundary. When normalized to the quantum of flux ^q, this 
electromagnetic flux is given by the following expression for the associated surface integral 
(the notation should not be confused with that for the electromagnetic scalar potential in 
S.I. units), 

^[C] = -7^ / dV^V- (43) 
2-K Js 

If the contour C and the surface S are purely spacelike, this quantity ^[C] does indeed 
measure, in units of ^q, the total magnetic flux through that surface. In the case of a purely 
timelike rectangular contour with fixed boundaries in space, the quantity <5[C] is in fact a 
measure of the electrostatic potential difference between those two points in space. The pur- 
pose of this general analysis is to make manifest once again the Lorentz covariance properties 
of the discussion, which will eventually be restricted to time independent configurations in a 
planar geometry only. 

In fact, by its very definition, the flux $[C] is not only a Lorentz scalar, but also a U(l) 
gauge invariant quantity. Indeed, through Stokes theorem, we also have 

^C] = -^<f^du^a^, (44) 

which is recognized to be precisely the U(l) gauge invariant Wilson loop associated to the 
contour C. Using now the general first-order relation in (|39|), we obtain 

Note that this relation, when considered for all possible contours C in spacetime, provides 
the integrated global expression of the covariant London equations ( |37| ) which in their local 
differential form cannot encode the topological information which we are about to characte- 
rize. Working in terms of this global integrated form of the covariant London equations is 
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crucial to the resolution of the LGH equations in a way which also accounts for the possible 
topological vortex structure of its solutions. The values $[C] for all possible contours C in 
spacetime thus also provide the "global boundary conditions" connecting the gauge invariant 
and gauge dependent variables and their semi-decoupled system of differential equations in 
the manner described in the previous section. 

Since the angular variable 9{t,u) need not be single- valued throughout spacetime, the 
contour integral of its spacetime gradient df^9 could a priori lead to an arbitrary real winding 
number L[C] in correspondence with the specific contour C being considered, such that 



1 

2ttJc 



— j> du>' d^e = -L[C]. (46) 



Hence for an arbitrary contour C, the integrated (covariant) London equations read, 

^C] = L[C] - ^£dn^j^ = L[C] + l-j^du^^l^J^. (47) 

In the case of a contour C which only lies in space directions, this general expression reduces 
to the following relation for the corresponding magnetic flux, 

^C] = L[C] + ± £ • J = L[C] -1-j^du-j^. (48) 

When considering a contour C at infinity in the case of an unbounded superconducting 
domain, this expression leads to the well known property of magnetic flux quantization at 
infinity, since the electromagnetic current density contribution then vanishes while L[C] must 
be an integer for the topological reason recalled in the Introduction. 

In fact for all choices of contours, including those at a finite distance, the values for 
L[C\ are restricted by the necessary single- valuedness of the order parameter ij), which thus 
constrains the lack of single-valuedness in the phase parameter in correlation with the fact 
that the function / is real but with a sign which could vary throughout spacetime, while its 
overall sign is left unspecified. As explained in the Introduction, to best consider that issue 
let us now imagine taking the scalar field ij) around the closed contour C. Given the winding 
number L[C\ associated to this contour as defined in (46), under transport around C the 
phase variable 9 is shifted by the quantity (— 27rL[C]), 

9^6- 27rL[C7] , e'^ ^ e-'^i^L[c] ^ie ^ ^^g^ 

In turn, since the order parameter if) — which is assumed not to be vanishing throughout 
spacetime — must be single-valued under such a transformation, the function / itself must 
then transform according to 

V'^V , /^e^'-^t^l/. (50) 

However, since the function / is to take real values only (but not necessarily positive ones 
only), consistency of such transformations thus restricts the possible winding numbers L[C] 
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for any contour C to integer as well as half-integer values only. In the case of integer winding 
number values, the function / recovers its original sign after transport around the correspon- 
ding closed contours, while the phase 9 is then shifted by an even multiple of vr. In the case of 
half-integer values however, / changes sign under such a transformation, in correlation with 
the shift by an odd multiple of vr in the phase 6. Nevertheless in either case, the complex 
scalar field ip remains single-valued throughout spacetime, as it should. Note also that if 
one artificially restricts the function / to take positive values only, as is done in the usual 
classification theorems Q, only integer winding number configurations survive the analysis 
with the possibility of a vanishing order parameter at discrete locations only but not in a 
continuous fashion such as inside some domain walls. 

Given this topological characterization of solutions, in the case of a half- integer winding 
number L[C] the function / must necessarily vanish and change sign an odd number of times 
when taken around the closed contour C, while in the case of an integer winding number / 
may vanish and change sign an even number of times, including zero of course. Moreover in 
the case of a non vanishing winding number L[C] — be it integer or half-integer — for a contour 
C which shrinks to a point, necessarily the order parameter ij) and thus the function / must 
vanish at that point, since the phase variable 6 is then ill-defined at that point being shifted 
by (— 27rL[C]) when taken around that point. Hence by continuity, given a contour C of non 
vanishing winding number L\C] and a surface S with C as a boundary, the general solutions 
to the LGH equations may be characterized by having the order parameter tp vanish at some 
point on S as well as on a series of continuous lines lying within S and emanating from that 
point, the number of such lines being even in the case of an integer winding number and 
odd for a half-integer one. Viewed in three space dimensions, such configurations correspond 
to vortex solutions extending along some spacelike one dimensional curve and to which an 
even or odd number of two dimensional finite or semi-infinite domain walls are attached at 
one of their edges according to whether the value of the vortex winding number is integer or 
half-integer, respectively^ 

In order to characterize the possible network of vortices and domain walls specific to a 
given configuration, let us turn again to the gauge invariance properties of the abelian U(l) 
Higgs model. As emphasized previously, the flux $[C] is both Lorentz and gauge invariant. 
In fact, under gauge transformations, see (]4^, the winding number contribution to $[C] 
transforms as, 

while the contribution from the current is gauge invariant by itself. Hence, it would appear 
that the winding numbers of any solution could be shifted away through some appropriate 
gauge transformation x whose contour integrals exactly cancel the winding number contri- 
bution to the flux ^[C] for all contours C. Such gauge transformations however, leading to 
topology change in the gauge and Higgs fields, are not allowed. Indeed, they would neces- 
sarily correspond to functions x which are not well defined throughout spacetime, possessing 

^Note that in the general case, the winding number L[C] is a function of time for any specific fixed contour 
C in space. This function is piece-wise constant, with integer or half-integer discontinuities at those specific 
instants when integer or half-integer vortices leave or enter the closed contour C. 
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singularities precisely canceling those of the vortex solutions. Since only gauge transforma- 
tions X which are regular and single-valued throughout spacetime are acceptable, the winding 
numbers L[C] as well as the fluxes $[C] are indeed gauge invariant physical observables for 
whatever choice of contour C. 

On the other hand, for solutions of non vanishing winding number, the phase variable 
must necessarily possess specific singularities at some points in spacetime, since some of 
its contour integrals are non vanishing even when shrinking to a point which then corres- 
ponds to the location of a vortex. Nevertheless, the regular component of 6 may be changed 
at will through arbitrary gauge transformations (^) parametrized by regular functions Xi 
without changing the winding numbers of but affecting nonetheless the gauge potential 
accordingly. This remark thus suggests that it may be possible to gauge away entirely any 
regular contribution to the phase variable 6 of the order parameter ij), leaving over only the 
singular contributions responsible for all winding numbers encoded in while modifying the 
gauge potential appropriately through the relevant gauge transformation x- To show how 
this is can be achieved, let us assume to have constructed a specific function whose winding 
numbers for all spacetime contours C reproduce exactly all those of 9 given a specific vortex 
configuration solving the LGH equations (such a function will be constructed explicitly in 
the next section). Since the difference (^o ~ ^) is necessarily a regular function which is 
well defined throughout spacetime, let us then apply the gauge transformation of parameter 
X = Oq — 0, leading to 

e' = eo , a'^ = a^-d^e + d^do. (52) 

However, given the second-order equation in (|3^), the transformed gauge potential a'^ is such 
that 

df'a'^ = d^a^' - d^d^^e + d^d^^Oo = d^f + a^a^^o- (53) 

Finally, since the current is gauge invariant, this last relation for a'^ remains invariant 
under all further gauge transformations x such that dfj_d^x = 0- Whether only the trivial 
solution X = obeys this equation depends on the choice of boundary conditions for the 
electromagnetic sector of the system at infinityP]. Whatever the case may be in this respect, 
the important conclusion is thus as follows. Given a specific configuration of bound vortices 
and domain walls whose topological structure is encoded into the function one may impose 
the "vortex gauge fixing condition"^, 

df.a'^ = d^f + 5^9^00- (54) 

The solution for the gauge dependent variables of the system is then given by 

= 6o , af, = df,eo + j^. (55) 

Depending on the choice of boundary conditions at infinity, the vortex gauge fixing ( ^4|) 
may not be complete, in which case both 6 and are defined only up to those gauge 

^°We refrain from being more specific on this point since one often considers the application of homogeneous 
external electric and magnetic fields outside the superconductor. 

^'^Note that this gauge fixing condition is Lorentz invariant, and that for time independent axially symmetric 
configurations, it reduces to the Coulomb-London gauge fixing condition d ■ a — (see the next section). 
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transformations whose parameter function x ^Iso obeys the equation 



d^d^x = 0, (56) 

but leading nevertheless to a construction which is consistent with the above expressions 
for the solutions for both 9 and a^, since both these quantities are then gauge transformed 
accordingly. Hence, the vortex gauge fixing condition ( [5^ ) together with the knowledge of the 
function 9q does indeed completely solve the LGH equations in the sector of gauge dependent 
variables, as announced previously. Note that when non trivial gauge transformations such 
that dfj^d^x = exist, the arbitrariness that such a situation leads to in terms of 6 and may 
in fact be absorbed into the choice of function which then also suffers the same physically 
irrelevant ambiguity. 



2.4 The function 6q and the double sheeted covering of the plane 

In order to understand how to construct the function in the general situation, let us 
first restrict to time independent configurations in the planar case. In other words, we 
assume that the system has been dimensionally reduced to two flat space dimensions, with the 
time independent fields dependent only on the two coordinates of that plane, the magnetic 
field b purely transverse to that plane and the electric field e and charge distribution 
vanishing. Choosing coordinates on the plane — normalized to the length scale A — which are 
either cartesian, with x = u^, y = u^, or polar, with u = y/ x'^ + y^, (j) = Arctan(y/x) and the 
specific evaluation — vr < 4> < +7r, it is useful to introduce the complex combinations 

z = x + iy = ue*"^ , z* = x - iy = ne"*"^. (57) 



Let us now consider a specific collection of K vortices described as follows. Each 
of these vortices, labelled by k = 1,2, . . . , K , has integer or half-integer winding number 
Lk = Nk/2, Nk being thus an even or odd integer, respectively. The total winding number 
of the configuration is thus L = N/2 with L = J2k=iLk and N = J2k=i-^k- In addition, 
each of these vortices has a position in the plane corresponding to the complex parameter 



Zk = Xk + iyk 



UkC 



l<Pk 



In order to construct the function 9o associated to this vortex 



configuration, let us also introduce for each value of k = 1,2,...,K a collection of 
complex functions 9k,nk{z) labelled hy Uk = 1,2, Nk and dependent only on the complex 
variable z. Then, a choice of function 9o{x,y) = 9o{u,(j)) which reproduces all the contour 
integrals associated to this vortex configuration is given by the following expression. 



1 K Nk 

k=l nir=l 



Zk) 



+i9t 



(58) 



where a specific evaluation of the complex logarithmic function is to be chosen (for example 
with its branch cut along the real negative axis in the plane, namely at </> = zbvr for ln(z)). 

The meaning of this construction is as follows. The introduction of the index taking 
a total of values is in fact related to the interpretation discussed in the Introduction 
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which views all possible configurations as built up from 1/2- vortices and 1/2-domain walls 
which are combined and stacked on top of one another. Each value for the index thus 
refers to each such 1/2- vortex as a basic building block, with in particular the corresponding 
function ^^.n^ i^) specifying the position in the plane (x, y) of the logarithmic branch cut in 
the quantum phase = Oq which is associated to that particular 1/2- vortex and whose branch 
point is at the position of that vortex. However, note that one may also write, up to some 
specific integer multiple of tt to be added on the r.h.s., 



1 



K 



In 
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Zk 



2; — zl 
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k,nk 



k,nk 



{Z* 



(59) 



Obviously, since the double sum on k and ra^ in the r.h.s. of this expression defines a real 
function x{x, y) which trivially satisfies the equation dfj,d^x = 0; this term may be gauged 
away altogether, leaving over only the first simple sum over k. In other words, the positions 
of the different logarithmic branch cuts in the function 9 = 9q may be moved around in the 
plane {x, y) through arbitrary gauge transformations; only their branch points are fixed at 
the 1/2- vortices locations. In particular, these branch cuts may be chosen to lie exactly along 
the lines of vanishing order parameter within the domain walls that may exist. Hence, when 
taking that degree of freedom into account which allows to gauge away all functions ^^^^^(z), 
there is no loss of generality in the following specific choice, for instance, for the function 
associated to such an arbitrary vortex configuration, 
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?o = XI ^k In 
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(60) 



As a particular example, consider a single vortex of winding number L at the center of 
the plane. One then has 

so that the order parameter takes the general form 

^{u,^) = f{u,(l>)e-''^^. (62) 

Having chosen to work with the angular range — tt < (p < +7r, clearly the phase factor 
9 = 60 = —Lcf) of this configuration has a branch cut at ^ = ibTr starting at the origin, 
while f{u,(p) must then change sign on that branch cut for a half- integer winding number 
L. A similar discussion applies in general, by considering the neighbourhood of any vortex 
in the plane. Hence, whenever half-integer vortices are present — which, as discussed in the 
Introduction, should indeed be the generic situation — , these branch cut properties show that 
the function / is in fact defined over a double covering of the plane, or a finite domain of 
it in the case of a bounded superconductor, with branch points at the positions of these 
vortices and branch cuts on those lines on which the order parameter vanishes within the 
domain walls which are bound to these vortices. The function f{u,(f)) then changes sign 
on these branch cuts, in such a way that it be continuous on the specific double sheeted 
covering of the plane which is associated to the considered configuration of bound vortices 
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and domain walls. Clearly, the same general picture also applies to integer vortices bound 
onto domain walls. It is only in the case of isolated integer vortices, thus not bound onto the 
edges of any domain walls — namely, the only situation thought to be possible until now — , 
that the function f{u, 4>) is in fact continuously single- valued on the plane (x, y) itself, rather 
than on some double covering of it, in which case it is indeed justified to assume that it 
be always positive (or negative) everywhere, with the exception of isolated points where it 
vanishes. Nonetheless, given the arbitrary overall sign for the function /(«,</)), it remains 
more appropriate, even in such a case, to still view that function as being defined on a double 
sheeted covering of the plane in which the branch cuts between the two sheets have 

degenerated in a continuous fashion into single points which then correspond to the positions 
of the isolated integer vortices. 

Hence by allowing all the degrees of freedom hidden in the polar parametrization of 
the order parameter, tp = /e*^, to manifest themselves in a way made consistent by the U(l) 
gauge invariance of the system, we see that the general solution to the LGH equations is 
associated to a specific double covering of the plane or a finite domain of it, in the manner 
just described. In particular, the unspecified overall sign of the function / is then also 
seen to correspond to a choice of sheet in this double sheeted covering of the plane. When 
only isolated integer vortices are involved, these two sheets become almost disconnected 
by touching only at isolated points, thereby allowing the sign of / to be fixed to remain 
either positive (or negative) as has always been assumed implicitly until now. But when 
domain walls are involved as well, then necessarily the two sheets become intertwined in a 
topologically non trivial manner precisely along these domain walls. Finally, in the case of 
a planar domain of finite extent, some of the branch cuts in / may extend all the way up 
to the boundary of the domain, and this general picture is compounded even further by the 
possibility of annular current flows surrouding the vortices, and thus adding further concentric 
closed branch cut structures in an almost periodic fashion to the double covering of the planar 
domain, with the function / changing sign at the boundaries of successive annular flows. 

As a last remark concerning the function as constructed above for time independent 
planar configurations, note that it satisfies the equation 5^9^00 = 0- Hence in such a situation 
the vortex gauge fixing condition (^^ reduces to the condition that df^a^ = dfj_j^, namely 
dxo} + dyO? = dxj^ + dyj'^. For an axially symmetric time independent configuration, we 
then have the Coulomb-London gauge fixing condition d ■ a = 0. 

Let us now return to the general situation of a collection of integer and half-integer 
vortices bound onto the edges of a collection of domain walls, all in a time dependent fashion 
and moving in a three dimensional space. In such a case, it is possible to define^ a curvi- 
linear system of coordinates associated to a specific foliation of spacetime, such that two of 
the coordinates define spacelike surfaces which are locally transverse to each of the vortices 
throughout space, with the third spacelike coordinate locally transverse to these surfaces, 
and finally the fourth timelike coordinate transerve — in terms of the Minkowski metric — to 
all three spacelike coordinates. In other words, at each instant in time it is possible to view 
the vortex system as being obtained from the previous planar description, by appropriately 
bending these planes, and by stacking them on top of one another in a transverse direction 

^^Such a coordinate system is in fact not unique. 
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in such a manner that all vortices and domain walls remain continuous throughout space. 
The coordinates in the planes locally transverse to all vortices then play the same role as 
the cartesian or polar coordinates considered in the planar case above. Hence, it is possible 
to construct again the function in this general case using exactly the same functional ex- 
pression as above, but this time in terms of the curvilinear coordinates associated to these 
transverse planes. The only difference is then that the complex coordinates -Zfe(T) and func- 
tions 9k,n^{z,T) specifying the positions of the vortices and of the logarithmic branch cuts of 
6 = Oq, become functions of time as well (with the functions Ok^n,,{z,T) amenable to being 
gauged away as before). This is thus how the topological structure and time dependency 
of a given vortex and domain wall configuration may be encoded into the function Oq^ from 
which the solutions for = Oq and = d^O^ + may be obtained in the vortex gauge 
dutt'^ = dfij^ + dfj,d^9Q in the manner described previously. Note that in this general case, 
the function 6q no longer obeys the equation d^d'^Oo = 0. 

This concludes the general discussion of the topological classification of all solutions to 
the coupled Maxwell and LGH equations. The sector of gauge dependent variables, namely 
9 and a^, is determined in terms of a function possibly defined up to some gauge trans- 
formation, which encodes the gauge invariant topological structure characterized through all 
the winding numbers which are associated to a given configuration of vortices and domain 
walls. The sector of gauge invariant variables, namely / — up to an overall sign — , and j, is 
determined by solving the coupled LGH and Maxwell equations ( |3^ ) and (33), subjected to 
specific boundary conditions which include the "global boundary conditions" obtained from 
the electromagnetic flux values $[C] for all possible closed contours C in spacetime, which in 
fact provide the global integrated London equations and depend on the vortex configuration 
through its winding numbers. This last set of equations and boundary conditions is totally 
gauge invariant, which is an advantage when solving this system. Furthermore, one should 
keep in mind that the function / is in fact defined on a double sheeted covering of the plane, 
with branch points at vortices and branch cuts along the surfaces of vanishing order param- 
eter lying within the domain walls. Finally, the (free) energy of such configurations may be 
obtained from the general expression (|40|). 

In this latter respect, it is interesting to consider a few sample situations corresponding 
to a single straight vortex of winding number L placed at the center of the infinite plane, 
for which we have seen that the winding phase factor is given by e-'^^. In the case of an 
isolated integer vortex, the physical configuration is in fact invariant under rotations around 
the vortex, since indeed only the order parameter 'ip(u,4>) = /(n)e~*^'^ is not cylindrically 
symmetric because of its winding phase factor, while no domain wall is attached to the vortex 
which would otherwise break this rotational invariance. Consequently, at a fixed distance u 
from the vortex, the order parameter retains a fixed vacuum expectation value |/(u)| < 1, 
so that when taking the order parameter around the contour of radius u from cp = — vr to 
(/) = +Tr, the values reached by ip in the complex plane all lie on the circle of radius |/(u)| < 1 
which is then wound L times starting from ip = e™^ f{u) = (— l)^/(n) and ending back at the 
same point on that circle. This circle also lies somewhere in between the absolute minimum 
IV'I = 1 and the local maximum = of the surface in the shape of the bottom of a wine 
bottle which is defined by the Higgs potential (1 — IV'P) ■ 
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Consider now again an integer vortex, say with L = 1, but this time bound to an even 
number of domain walls. Rotational invariance being then broken, when going around the 
contour of radius u the order parameter tl^{u, cj)) = f{u, 0)6^*^"^ no longer retains a constant 
norm. Rather, starting from one of the domains walls where ip = 0, f = 0, the order 
parameter moves in a continuous fashion from the top of the Higgs potential down to some 
maximal value for the norm IV'I < 1 and back to the local maximum at = by following a 
closed path which passes through the latter point as often as there arc domain walls bound to 
the vortex. More specifically, when there are two domain walls, a single such loop is followed 
from top to bottom and back twice in the same direction. But when there are four domain 
walls for instance, then the entire closed circuit lying on the surface of the Higgs potential 
has the shape of a four-leaved clover whose center is at the top of the Higgs potential, and 
so on for still larger even numbers of domain walls. 

Similarly in the case of L = 1/2 and a single domain wall, the order parameter ip(u, (p) = 
f{u, (j))e-'^'t' follows only once a closed loop running from the top to some lowest point on the 
Higgs potential surface and back. When there are three domain walls bound to the L = 1/2 
vortex, a three-leaved clover shape is obtained, and so on. Taking a L = 3/2 vortex to which 
three domain walls are bound, again a single loop running from top to bottom and back is 
followed three times, while when only one domain wall is bound to the L = 3/2 vortex the 
closed circuit is followed only once and reaches the top of the Higgs potential only once, but 
it then also explores surroundings lying on opposite sides of the top of the Higgs potential by 
following a spiraling path which leaves and returns at the top in a direction perpendicular to 
that of its lowest point. 

In fact, for symmetry and energy reasons, one would expect that for half-integer vortices 
of winding number L = N/2 the configuration of lowest energy would cither be that with all 
|A^| bound domain walls stacked on top of one another, or else spread out in a maximally 
symmetric star-like shape with an angular opening of 27r/| between each. For a given vortex 
of integer winding number L = N/2, there also exists the third possibility that no domain 
wall whatsoever is bound to the vortex. Presumably, which of these different possibilities 
is actually realized in each of these cases would also depend on the value of the scalar self- 
coupling K or Aq. Indeed, it is for either one of these numbers of domain walls that the 
order parameter, when taken around a closed contour surrounding the vortex, would explore 
the more regions closer to the bottom of the Higgs potential, thus leading to a smaller 
condensation energy contribution. Only a detailed numerical analysis however, would resolve 
this issue. 

3 Two Dimensions and BPS Bounds 
3.1 Two dimensional reduction 

For the remainder of the paper, we shall restrict the discussion to time independent configu- 
rations in the plane, in the manner already described previously. Furthermore, since when 

considering bounded domains in the plane these will have the topology of cither a disk or an 
annulus, the appropriate choice of coordinates is the polar one with the u and — tt < (p < +7r 
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variables also introduced previously. In particular, a disk topology will have a radius b, 
namely Uf^ = b/X given our choice of normalized units, while an annular topology will have 
inner and outer (normalized) radii Ua = a/A and Ub > Ua, respectively. 

An arbitrary time independent configuration in the plane is represented in terms of 
the following set of quantities. The electric field e, scalar gauge potential if and charge 
distribution all vanish. The magnetic field b is perpendicular to the plane with a sin- 
gle component b{u,(j)), such that h{u,4>) = (0, 0, 6(ti, 0)). Finally, the vectors j{u,(j)) and 
d{u,(p) only have radial and azimuthal components, namely j{u,(j)) = {ju{u,<j)), j(j,{u,(j)),0) 
and a{u,<j)) = {au{u,(j)),a(j,{u,<j)),0). All these non vanishing quantities are functions of u 
and (j) only, but neither of the coordinate transverse to the plane nor of time, as are then 
also the order parameter ijj{u,(j)) and its polar decomposition variables f{u,(j)) and 9{u,(j)). 
In particular, note that this restriction on the decomposition of the vector potential a is 
consistent with the vortex gauge fixing condition (|5j). As a matter of fact, it also proves 
useful to introduce the following notation, 

j{u,(l)) = ju{u,(p) , g{u,(l)) = uj^{u,(p), (63) 

since these quantities appear naturally in all expressions. 

Given such specific symmetry restrictions — namely translation invariance in time and 
in the direction transverse to the plane, or in other words parallel to the straight vortices and 
domain walls — , it is immediate to determine the form of the ensuing system of equations. 
For the LGH equation (52|), one finds 



nduudu + dl]f= [{ujf + g^]f- k^u^ (l - f) f, (64) 
while the inhomogeneous Maxwell equations ( |33[) reduce to (in first-order form) 

uduh = fg , d^b = -fuj, (65) 
where the magnetic field (transverse component) is given by the local London equation 

b = -dug - -d^j. (66) 
u u 



The current conservation equation, which follows from (|65| ) of course, reads, 

udu (fuj) + [fg) = 0. (67) 

The set of coupled differential equations (64) and (pBf) is subjected to a series of boun- 



dary conditions. First, the "global boundary conditions", corresponding to the global London 
equations, are 

^C] = L[C] + ^ i du-j, (68) 
ivr Jc 

where C is any possible choice of closed contour C in the plane, L[C] is the winding number 
of the quantum phase variable 6{u,(j)) = 9o{u,(l)) (in the vortex gauge) associated to that 
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contour, and $[C] is the magnetic flux — normahzed to the quantum of flux $o — through the 
surface S in the plane bounded by C, 



^C] = — I dud(t>ub{u,^). (69) 
2iT Js 

In particular, for a circular contour centered on the origin n = of the plane and of radius 
uq, we thus have the condition 

— / du 



2tt 



/ du d4>ub{u,(l)) = ^[uo]= L[uo] + — dcp g{uo,4>), (70) 

Jo J-n ^Vr J-7r 



where L[uo] stands of course for the associated winding number in 6{uo,(p). 

In the case of the infinite plane, the local boundary conditions at infinity which complete 
those in ( |68[ ) and which are required for finite energy configurations are such that /(n, </>) 
approaches either one of the two values / = ±1, while the quantities j{u, <j)) and g{u, (p) must 
both vanish, since they directly define the components of the electromagnetic current density 
J = —f^j in the plane. 

In the case of a bounded domain in the plane, there may exist an applied external 
magnetic field 6ext solely transverse to the plane with a component 6ext- Based on Maxwell's 
equations in vacuum outside the bounded domain, in the specific instance of such a geometry 
for 6(n, (j)) it then follows that the magnetic field b{u, (/)) must indeed retain a constant homo- 
geneous value throughout space in that region. Consequently, given the remaining boundary 
conditions relevant to f^j and df discussed previously, it then follows that in the case of the 



disk topology the required local boundary conditions which complete those in (65) are 



b{ub,(j)) = be^t , j{ub,(p) = , {duf) {ub,(j)) = 0. (71) 



In the case of the annular topology, the same boundary conditions apply of course at 
the outer boundary at u = Ub- Similarly at the inner boundary u = Ua, one must have 

b{Ua,^)=ba , j{Ua,(p)=0 , {duf){Ua,<P)=0, (72) 

where ba represents the value of the constant magnetic field (transverse component) within 
the hole of the annulus. The value for ba is to be determined from the total magnetic flux 
^[ua] through that hole which is to be expressed in terms of the angular integral of g{ua, 4>) 
and the winding number L[ua\ in the quantum phase 9[ua, 4>) of the order parameter 'ip{ua, 4>) 
on the inner boundary, as given in ( |70| ) with uq = Ua- Note also that the fact that b{u, (j)) must 
be (^^-independent on the boundaries is consistent with the determination of b{u, (j)) inside the 
superconductor in terms of j{u,(j)) and g{u,(j)) and with the boundary conditions on j{u,(l)) 
(see (II)). 

Finally let us consider the expression for the free energy of the system. Due to the 
invariance of the considered configurations under translations transverse to the plane, the 
three dimensional integral in (^) is in fact infinite. We should thus rather consider the 
expression for the free energy per unit of length in the transverse direction, while we may 
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also absorb the overall normalization factor in (^0|) in that choice, thus leading in those units 
to the expression 



£ = l^du £ dcpuf^ib- 6e.t]' + (duff + ^ {d^ff + f (i' + ^5') + (1 - • 

(73) 

Note that the constant term (— k^/2) has not been added to the integrand of this expression, 
keeping in mind the possible infinite planar topology in some applications (for the finite disk 
and annular topologies, this term will be included later on). Note also that even though the 
radial integration extends throughout the infinite plane, it is only within the domain of the 
superconductor that /, j and g do not vanish, while in the outer region of the supercondutor 
(if Ub is finite) we always have b = 6cxti so that in fact there is no contribution to the free 
energy for u > Ub in the cases of the disk and annular topologies (there is one for u < Ua in the 
annulus case). Actually, when evaluated for a specific solution to the differential equations 
( |64[ ) and (|65|), the expression for the free energy as given in (73) also reduces to 



poo pn f 11^ 

The above set of equations and boundary conditions thus defines the problem to be 
solved. The only information missing is that which provides the topological structure of the 
vortex and domain wall configuration which is being considered. This structure is characteri- 
zed by the set of winding numbers L[C] which appear in the global boundary conditions (|6^). 
The gauge dependent details of that topological structure are specified through the choice of 
function Oq{u,4>)^ which in the vortex gauge ( [5^ ) may be expressed as in ( pO| ) as explained 
previously, and in terms of which the quantities 9{u,4>) = 9o{u,4>) and a{u,<j)) = j — dOo are 
then obtained. Nevertheless, it is only the gauge invariant content of ^o, through the winding 
numbers L[C], which is involved in the resolution of the remaining quantities f{u,(l)), j{u,4>) 
and g{u, <j)), which are to be determined from the above differential equations and boundary 
conditions, while also keeping in mind the double sheeted covering of the plane which is 
associated to the vortex and domain wall configuration being considered and over which the 
function /(u, </>) is to be defined. 

As a final point, consider the case of the infinite plane or the disk with a vortex of 
winding number Lq 7^ placed exactly at its center u = 0, and a circular contour of radius uq 
surrouding that vortex. In the limit that the contour shrinks to the origin, since the magnetic 
fiux ^[uq] must vanish — the magnetic field being finite at u = — , the azimuthal component 
of the current density J = —f'^j at that point must be such that 

/ (i(t>9(.uo,4>) = -Lo, (75) 



which necessarily requires 



5(0,0) = -Lo. (76) 



On the one hand, since j(j){u, (f) = g{u, 4>)/u, this result shows that the current j is necessarily 
ill-defined at the location of vortices. But on the other hand, since the order parameter ip 
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and the function / vanish at the location of the vortex, while the azimuthal component of the 
electromagnetic current J must vanish at the same point (this is a consequence of the polar 
coordinate parametrization which is being used), this result also shows that the singularity 
in j = — J//^ is just mild enough to be screened by the zero in il) in such a way that the 
physical current J remains well-defined and finite throughout space as it should, even at the 
location of vortices. 



3.2 BPS bounds 



The resolution of the LGH equations, even in the plane, requires a numerical analysis, since no 
analytic solutions are known. Nevertheless, it is possible to establish some general properties 
for the possible solutions in some particular situations. Namely, for specific values of the 
scalar self-coupling k or Aq, it is possible to establish BPS lower bounds |^ on the energy 
of vortex configurations, as recalled in the Introduction. 

To this aim, let us consider the free energy expressed in ( [73| ) in the case that no external 
magnetic field is applied, 6ext = 0, and for the disk and infinite planar topologies^. Through 
integration by parts and the identification of the geometric meaning of the induced surface 
terms, one then finds that 



£ 
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Lk + §9ndu-{l-nj 
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duf + v^fa 



vfj 



(77) 

where rj = ±1 is some choice of sign to be specified presently, il. stands for the planar domain 
in which f ^ 0, namely that of the superconductor, and d^l for the boundary of that domain. 



Finally, L 
the domain 



Lk is the total winding number of the vortex configuration lying within 



The remarkable feature of this relation is that, with the exception of the surface term 
at do, and the very last term which depends on the self-coupling k, this expression gives 
the energy of any solution in terms of a sum of integrated positive quantities which only 
involve first-order variations — as does also the expression in ( ^3|) — but in such a way that 
the topological content of the configuration is also made explicit through the total winding 
number L. Furthermore, when attempting to minimize the value of the energy by setting 
to zero the positive integrated quantities, thus leading to the following first-order differential 
equations, 



b = -dug - -d^j = ^r](l- f) , duf 
u u 2 \ ' 



1 



1 



d4>f = Vufj 



(78) 



it is immediate to verify that these equations then also imply the second-order equations in 



(p4D and (pq) , but only if the coupling takes the specific value k = Kc with Kc = 1/ v 2. In other 



The annulus case is left to the reader. 
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words, when the scalar self-couphng takes the critical value k = k^, the system of first-order 



equations (78) integrates the second-order LGH equations of motion (the issue of boundary 
conditions will be addressed shortly). Note also that it is for the same critical coupling Kc 
that the very last term in (^) does not contribute to the energy of such configurations. 

In order to fix the choice of sign for ry = ±1, let us note that the lowest possible value 
taken by £ as defined in (|7^ ) is zero, corresponding to the trivial solution b{u, (j)) = 6ext (here 
^ext = 0), f{u,(l)) = ±1, j{u,(j)) = and g{u,(j)) = in the absence of any vortex in the 
plane. Consequenly, all other solutions must necessarily possess a strictly positive energy 
value, which dictates the following choice of sign 

r? = sign VLfc =signL. (79) 



Hence, given this choice for r/ and the specific value k = Kc, the energy of any solution to the 
first-order equations (^) — which then also obey the second-order ones — is given by 



£ = 2Tr 



K 

k=l 



+ {s\gnL) j>^Ju-(l-f)l (80) 



However, it is only in the case of the infinite plane that the surface term on the boundary 
dil does not contribute, since the boundary conditions are then such that both (1 — /^) 
and j vanish at infinity. Furthermore, it is clear that these boundary conditions are also 
consistent with the first-order equations (|78|). Hence finally, one concludes that all solutions 
to the first-order equations ([78|), which also obey the second-order LGH equations and the 
appropriate boundary conditions in the infinite plane when k = Kc, saturate the BPS value 
which determines their energy or mass in terms of a topological invariant, namely the total 
winding number of that configuration. 



£ = 2tt 



K 



k=l 



27r|L|. (81) 



Conversely, it has been shownQ that under both the assumptions of an everywhere positive 
function / and possessing only a discrete set of zeroes of integer positive degree, all the 
solutions to the second-order LGH equations and boundary conditions in the infinite plane 
with K, = Kc also solve the first-order equations ([T^) with t] = signL. In other words, for the 
critical scalar self-coupling k = Kc, any solution to the LGH equations in the infinite plane 
with isolated integer vortices only and of fixed total winding number L saturates the BPS 
value above, thus also showing|jl5| that such vortices do not possess an interaction energy 
which would depend on their relative positions. 

However, when relaxing these restrictions on the function /, and thereby accounting 
for the possibility of half-integer vortices and bound domain walls in the infinite plane, one 
must conclude that such configurations cannot solve the first-order equations (|78|) and thus 
cannot saturate the BPS value (]8l|). Indeed, the first of these first-order equations would 
imply that on the surface of vanishing order parameter, / = 0, within any domain wall, the 
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magnetic field b would retain a constant value of t]/2 irrespective of the distance between the 
vortices which are bound onto the edges of that domain wall and irrespective of the winding 
numbers of these vortices. Such a property seems very unlickely, not only because of this lack 
of dependency of the magnetic field on the separation between the vortices at the edges of 
domain walls and on their winding numbers, but also because the magnetic energy density 
contribution to £ would then grow linearly with that distance, a fact which would clearly be 
inconsistent with the saturation of the BPS value. On the other hand, we also know that 
the further a domain wall is streched (provided the vortices bound onto its two edges no 
longer overlap), the more the condensation energy stored into the domain wall contributes 
to £ in an almost linear fashion as well, again a property which would be incompatible with 
the saturation of the BPS value (^l]). Hence, based on such physics arguments, it must be 
concluded that even for the critical value k = Kc, only isolated integer vortex solutions to the 



LGH equations in the infinite plane both solve the first-order equations (78) and saturate the 
BPS value, while any other vortex configuration including thus domain walls cannot share 
these properties. Nonetheless, such configurations with domain walls bound to vortices in 
the infinite plane do obey the following BPS strict lower bound 



£>27r 



K 

J2Lk =27r|L|. (82) 

k=l 



Indeed considering again ( |77D , the surface term at infinity vanishes in the infinite plane while 
the first two-dimensional volume integral is strictly positive since such configurations cannot 
obey the first-order equations ( [7^ ) (the last term vanishes for k = Kc). 

Let us now turn to the disk of finite radius u^y. In this case, one is led to the conclusion 
that the BPS lower bound £ > 27r|L| does not apply for solutions to the second-order LGH 
equations with k = Kc, even for isolated integer vortices. The first-order equations ( [78|) 
are in fact incompatible with the required boundary conditions (|7l[) |p!8[| . Indeed, even when 
ftext = 0, these equations and boundary conditions together imply that 

/(ub,0) = ±l , g{ubA) = ^ , {d^f){uf„cP) = 0. (83) 

Clearly if J2k=i = L ^ 0, such conditions can be met only if the boundary of the disk is 
at infinity, Ub oo, thereby recovering the previous discussion in the infinite plane. This 
conclusion on its own however, would still not exclude the relevance of the BPS lower bound 
to solutions in the disk. In fact, even had the first-order equations been consistent with the 
boundary conditions, there cannot exist a BPS lower bound on the energy of any solution in 
a finite disk, including the case of isolated integer vortices. Indeed, the surface contribution 
du ■ {1 — f'^)j does not vanish in such a situation, and could a priori carry whatever sign 
depending on the details of the distribution of vortices within the disk. 

The situation with respect to possible BPS lower bounds having been understood in 
the case of the critical coupling k = Kc, let us now consider again the general expression { fl7\j 
relevant whatever the value for n. Clearly, for couplings larger than the critical one, k > Hc, 
and configurations in the infinite plane thus leading to a vanishing surface term at infinity, 
the energy always obeys the BPS strict lower bound in (^2|) whatever the solution to the 
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LGH equations. Indeed, on the one hand, the very last contribution in ( [77| ) proportional to 
(k^ — 1/2) is then always strictly positive while, on the other hand, the remaining integrated 
positive terms cannot vanish either since otherwise this would require the first-order equations 
( [78|) (with 1] = signL) to be obeyed by solutions to the LGH equations, which is excluded 
when K ^ Kc- In the case of the finite disk with k > Kc, again no BPS lower bound may be 
given because of the non vanishing contribution of either sign from the surface term on the 
disk boundary. Values for the energy £ for configurations in the disk with k > Kc could a 
priori be smaller as well as larger than the BPS value 2tt J2k=i = 27r|L| depending on 
the distribution of vortices. 

Finally, when k < Kc, the same conclusion as to the absence of a BPS bound must be 
drawn, whether for vortex configurations in the infinite plane or the disk, since in that case, 
the contribution of the very last term in ([77| ) is always strictly negative. 

3.3 Singular BPS states 

There exists however, a specific situation, albeit a singular one, for which all solutions to the 
LGH equations, including domains walls, do in fact saturate the BPS lower bound (^T|). As 
the above discussion has shown, one crucial property for saturating the BPS bound is that 
the last integral in (|77| ) vanishes identically |lO|. This may happen in either one of two 
ways, namely if k = Kc as already discussed, or else if = 1 throughout the plane except 
possibly on a subset of zero measure. Furthermore, if the latter situation is realized, then 
the energy becomes independent of k, including even an infinite value for that coupling. 

In the context of superconductivity, having chosen to normalize distance measurements 
to the penetration length A, the limit k — > +00 corresponds to a situation such that, compared 
to A, the coherence length ^ becomes infinitely small. Likewise, had we chosen to normalize 
length scales to the same limit would correspond to a situation such that, compared to ^, 
the magnetic field penetration length becomes infinite. Hence, given that the external field 
6ext is assumed to vanish in our present considerations, one should expect that in the limit 
K — > +00 the magnetic field b{u, (p) must vanish within the superconductor except possibly 
at a few discrete points. As we shall see, this is indeed what happens. 

In a particle physics context however, an infinite value for k would likewise imply either 
a vanishing gauge boson mass or else an infinite Higgs mass M/j. In the first instance, 
the vacuum expectation value a would have to vanish while sending Aq to infinity in order 
to keep Mfi non vanishing and finite, while in the second instance the parameter a is kept 
non vanishing and finite while Aq is also sent to infinity. Thus in either case, this limit 
corresponds to an infinitely strongly interacting Higgs sector, a regime which is incompatible 
with perturbative unitarity in the quantized field theory. Nevertheless, in the exploration 
of classical solutions to the equations of motion of the abelian U(l) Higgs model, the limit 
K +00 may still provide some useful insight. 

With this motivation in mind, let us consider then the second-order LGH equations in 
(|64[) and (p5|) subjected to the relevant local and global boundary conditions. In particular. 



since only the LGH equation (64) involves the parameter k, the limit n — > -|-oo may be taken 
in a consistent manner by first dividing that equation by k^. Consequently, the only allowed 
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solutions must be such that either / = Oor/ = iblat every point of the plane, precisely 
allowing for all types of vortex and domain wall configurations defined over a double covering 
of the plane or a finite domain of it. Considering then the remaining equations for b{u, (p) as 
well as the global boundary conditions (|68|), it then follows that given a vortex distribution 
of winding numbers and positions Uk {k = 1,2, ■ ■ ■ , K) in the plane or a finite domain of 
it, the solution to the LGH equations in the limit k +00 is given by 

K 

biu) = 27rY,L,6^^Hu-Uk) , jiu) = -dxbiu) , i; (u) = f (u) e''<^^^\ (84) 

k=l 

where the function f{u) thus takes the value / = at the position of vortices and domain 
walls, and otherwise either one of the values / = ±1 according to the double sheeted covering 
of the plane which is associated to that bound vortex and domain wall configuration. 

Obviously, such solutions are singular precisely at the positions of vortices and domain 
walls, with the magnetic field vanishing everywhere except at the vortices where it diverges, 
and the function / constant everywhere — on the double sheeted covering of the plane — except 
at its branch points and on its branch cuts where it vanishes discontinuously — thus displaying 
the singularities of the double sheeted covering of the plane. Furthermore, the current j and 
thus also the electromagnetic current J = —f'^j, vanishes everywhere, except again at the 
vortices where it is singular. The physical picture however, is clear. In the limit n — > +00, 
the coherence length ^ being zero or the penetration length A being infinite, the thickness of 
domain walls has shrunck to nothing while likewise vortices have been squeezed into single 
curves without any transverse extension but nevertheless each carrying the required finite 
value of magnetic flux, by having "sucked in" all the currents responsible for that magnetic 
flux. In other words, in the limit k — > +00, the system of vortices and bound domain walls 
has reduced to a collection, on the one hand, of membranes of zero thickness which carry no 
magnetic flux, are bound at their edges to Aharonov-Bohm flux lines of integer or half-integer 
flux or extend to the boundary in the case of a finite domain in the plane, and which, on 
the other hand, may possibly also appear in combination with some isolated Aharonov-Bohm 
flux lines of integer flux. Hence once again, such conflgurations are very much reminiscent of 



D-branes[22] and strings in M-theory, which are precisely such objects without any transverse 
extension. 



Let us now consider the expression (77) for the energy of these configurations. Since the 
magnetic field b, the current j and the Higgs potential (1 — IV'P)'^ vanish identically except 
on a subset of the plane of zero measure, all such configurations of Aharonov-Bohm flux 
lines and membranes of flxed total winding number L = J2^t=i indeed saturate the BPS 
lower bound, £ = 27r|L|, thus showing that all such solutions are indeed BPS states, albeit 
singular ones. In particular, this result also implies that such Aharonov-Bohm flux lines 
and membranes have no interactions with one another, whatever their relative positions. 
From the physics point of view again, this fact is rather obvious. The order parameter 
ijj reaches its constant vacuum expectation value {tpl = 1 and the magnetic fleld vanishes 
everywhere in space except at the positions of flux lines and membranes (the latter applies 
to {ipl only), so that indeed both the magnetic and condensation energy densities vanish 
identically irrespective of the distribution of flux lines and membranes in space. Note that 
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these conclusions as to the BPS character of the solutions to the LGH equations in the limit 
K — > +00 remain valid whether in the infinite plane or a disk of finite radius u^. The only 
restriction in the case of the infinite plane is that the total winding number L = J2^=i 
must then be an integer, without any membrane extending up to infinity. 

4 Half-Integer Vortex Solutions 

Having argued that the basic entities from which to build general solutions to the LGH 
equations are 1/2- vortices, 1/2-domain walls and annular current flows, in this section we 
consider the situation of a single vortex at the center of either a disk or an annulus of finite 
radius. Even though no exact analytical solutions exist, we shall try nevertheless to gain 
some insight into the nature of such solutions by using some approximations, and then turn 
to the results of a modest first attempt at a numerical analysis of half- integer vortices. 

When solutions with half-integer winding number are being considered, note that we 
shall only consider here the case of a single domain wall bound on such a vortex, and further 
assume that this domain wall lies along one of the radii of the disk or annulus. Indeed, even 
though such configurations necessarily break the axial symmetry of the disk or annulus, this 
breaking is kept to a minimum under such a restriction, which then remains covariant under 
rotations around the center of the disk or annulus which map such solutions into one another 
without changing their energy. This restriction is also in keeping with the fact that the 
natural tension which domain walls possess is such that in their lowest energy configuration 
domain walls must be straight not only in the direction transverse to the plane, but also 
within that plane. 

4.1 Annular vortices 

In order to argue for the possibility of annular vortices of integer as well as half-integer 
winding number, we shall use the same approach as that of Ref.[^ which uncovered the 
existence of such configurations in the isolated integer case by considering the emergence 
of these solutions from the trivial solution at the normal-superconducting phase transition. 
Namely, we shall restrict to the disk topology with a single vortex of winding number L at 
its center, in the absence of any external magnetic field. Clearly, the solution to the LGH 
equations which describes the system at the phase transition is then given by 

g{u,(f>) = -L , j{u,(t>)=0 , b{u,(t>) = , f{u,(t>) = 0, (85) 

so that the order parameter vanishes identically, 'ijj{u,(j)) = 0, with a vanishing energy for 
the state, £ = (for the remainder of the paper, the subtraction constant (— k;^/2) is being 
included in the definition of the free energy (|7^)), and a non vanishing winding number 
nonetheless. 

This solution may also be obtained by considering a rescaling by a real parameter /o 
of the order parameter, such that 

i^iuA) = h^{u,(t)) , fiu,(l)) = fof{u,<))), (86) 
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and then taking the hmit /o = in the LGH equations and relevant boundary conditions. 
Doing so, one finds that the function f{u, (p) satisfies the linearized LGH equation, 



uduudu + / = [{ujy +g'\f- K'u'f. (87) 

Given the solution in (|85|), as well as the requirement that f{u,(j)) should vanish at the origin 
when L ^ since this is the position of the vortex where f{u, cj)) must vanish, the linearized 
LGH equation possesses a single solution. The case of an isolated vortex of integer winding 



number has been discussed in Ref.|18]. Let us restrict here to that of a half-integer value for 
the winding number L of a single vortex which is bound onto a single domain wall lying along 
a radius of the disk (the generalization to more domain walls and integer winding numbers 
is immediate, as the reader will realize). 

Under these specific assumptions, and the fact that the equation for /(n, (p) in the limit 
/o = is linear, it is natural to consider the following separation of variables, 

/>,0) =/(n)sin(0/2), (88) 

which clearly displays a cut at i;^> = zbvr, having assumed the domain wall to be at (/> = 
where the branch cut for /(n, (p) in the double covering of the disk thus lies. The solution for 
f{u) is then simply a Bessel function of the first kind with index a = y/ L'^ + 1/4, 

f{u) = (^^y r{l + a) J J Ku), (89) 

where the normalization is chosen such that the lowest order term in a series expansion in 
u is f{u) ~ u". As is well known, such a Bessel function has an almost periodic oscillatory 
behaviour with an amplitude which asymptotically decreases as at large radii. 

Consequently, as soon as the parameter /o is turned on slightly away from zero, this 
oscillatory pattern hidden in the solution f{u,(j)) = will emerge from the vanishing con- 
densate, and lead to a pattern of concentric annular current flows of approximately constant 
width tt/k (in units of A) surrounding the half- integer vortex of winding number L which 
is bound onto the edge of a single domain wall (only the index a of the Bessel function as 
well as the periodicity of the trigonometric function in (p change when more domain walls 
are involved). Of course, the boundary conditions at u = Ub are then no longer satisfied, so 
that both the parameter /q as well as the remaining functions j and g must then be adjusted 
in order to obtain a solution of negative energy to the LGH equations, thereby also pushing 
outwards the pattern of concentric annular current flows without essentially changing their 



width PI 



Nevertheless, the existence of the Bessel function solution to the linearized LGH equa- 
tion at the phase transition provides a very strong argument for the existence in finite domains 
of annular vortex solutions of arbitrary half-integer winding number bound onto the edges 
of an arbitrary number of domain walls. The existence of such solutions without domain 



walls has indeed been established along such lines for integer vortices in Ref.[18|. There is 
no reason to doubt that the same conclusion should not also hold for integer and half-integer 
vortices bound onto domain walls. 
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Even though it is impossible to obtain an exact solution to the equations when turning 
on the parameter /o, even in the form of a power series expansion, it is interesting to consider 
the first-order corrections that this implies for the quantities j, g and b, without including 
the corrections implied for /(u, </)) itself. To lowest order in u, and still assuming only a single 
domain wall, one finds, when also imposing the current conservation equation, 

UJ{U,(p) - 2asin(0/2)' 

3K</>) = -L+l6on2-,^/oV("+^)sin2(0/2)-^-^ln(n/^o), (gg) 
b{u,cl>) = 6o-^/oV-sin2(,/./2), 
f{u,cf>) = fou^sm{4>/2), 

where bo stands for the value of the magnetic field at the vortex n = and uq for an arbitrary 
integration constant. These expressions are at least indicative of the singularities in j(u,(j)) 
and g{u, (j)) not only at the position of the vortex — a fact which was already demonstrated 
previously — but also on the surface of vanishing order parameter at </> = 0, namely on the 
branch cut in the function f{u,(p), but in such a way that the physical electromagnetic 
current J = —/"^j remains well defined and regular on this branch cut nonetheless, including 
the quantum tunnel effect through the domain wall in the proximity of the vortex as was 
described in the Introduction. These properties, announced previously, are thus made explicit 
in this approximate representation of a single vortex in the disk bound onto a single domain 
wall. 



4.2 The thin annulus hmit 



Another instance for which it is possible to gain insight into the nature of half-integer vortex 
solutions in a particular limit, is that when the width of the annular topology vanishes, the 
outer radius being kept fixed, namely the limit Ua ^ u^. Let us assume that a single winding 
number L is trapped in the center of the annulus of width Au = Ub — Ua and outer radius Ub- 
Consequently, the value ba of the magnetic field inside the hole must be such that 

— J d(t)g{ua,(l)) = -ulba- L. (91) 

In the limit that the annulus becomes infinitely thin, Ua — > Ub, the magnetic field 
both within the superconductor as well as inside the hole takes the value bext of the applied 
field. Moreover in that limit we shall also take the approximation that j{u,<j)) = 0, the 
annulus having no radial extent anymore, even though such a restriction may prove to be 
inconsistent with the current conservation condition. In this limit, not only does the above 
global boundary condition imply a specific restriction on the function g{ub,(l>) in terms of L 
and bext, but also the free energy ( [73| ) then reduces to the expression, 

£^^J[ # {(5<a/)' + fg' + I'^'nlf - K^lf^ . (92) 

In order to address the problem of minimizing the free energy and thus solve the LGH 
equations in the thin annulus limit, let first consider the case of an integer winding number 
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L without any domain wall lying along some radius within the annulus. Because of the axial 
symmetry of such a configuration, no (/)-dependency arises for the functions b, g and /, and 
one is left only with the following relations (with all these quantities evaluated at u = Ub, of 
course) 



9 



L , £: ~ 27r- 
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(93) 



It is immediate to show that the free energy is minimized for the following values 
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(94) 



In turns out that this quartic dependency of £ on the applied field 6ext provides a rather good 
approximation to exact numerical solutions for annuli of finite width Am small compared to 
Ub, since it presents a turn over point precisely at the values of 6ext where this expression 
vanishes. This property is very close to the behaviour of the exact solutions, while it is 
not reproduced by the quadratic dependency which is usually discussed in textbooks [|^, 17 1. 
Whatever the relevance of this latter remark, the important point is that when considered 
for all possible integer values L, the graphs for the free energy £ as a function of &cxt all cross 
one another in succession before their turn over points at £ = provided that nub > 1/2 and 
for values such that 
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k being an arbitrary positive, negative or zero integer. This, of course, is the experimentally 
observed behaviour as well, namely the Little-Parks effect ||l^, |l^, the crossing points occur- 
ring for a magnetic flux through the annulus equal to a half-integer multiple of the quantum 
of flux $0 (note also the lower bound on the annulus radius, b > ^/2). Furthermore, the 
absolute minimimum of the free energy is reached for 
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while the absolute maxima are reached for 
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(96) 



(97) 



Let us now consider a similar analysis in the case of a half-integer winding number 
L, which should thus also include some domain wall structure inside the thin annulus. The 
axial symmetry of the equations then being broken, all quantities b{<j)), g{(p) and /(</>) become 
(/)-dependent. Hence, the minimization of the free energy leads in this case to the differential 
equation, 

^f = fg^-^\l(l-f)f, (98) 
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while the function g[(p) is constrained by the condition 

^ £^ #5(0) = luX^t - L. (99) 

Finahy, let us also include the current conservation equation, which in the present limit 
reduces to 

d^fg) = 0. (100) 

Hence, we must have 

«W = 7^ • i|^ = ^-«'"'<^-^''^- 

with C some integration constant. 

First in the case that this constant vanishes, C = 0, we necessarily have g{(f)) = 0, 
which thus implies that such a situation is obtained only for values of the applied field such 
that 

-ulbe^t = k + - = L, (102) 

k being an arbitrary integer. These are precisely the values for which the energy of the thin 
annulus is degenerate for two consecutive integer winding numbers (provided nub > 1/2). 
Furthermore, the solution for the order parameter is then 

/(0) = ±tanh(^^<A), (103) 

leading to a value for the free energy of this configuration given by 
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(104) 



This value always lies above that £{k) for the integer winding number case at the same 
magnetic field values, whatever the value for KUb > 1/2. 

The solution found for /(</>) does indeed display a cut at (p = ibvr, in accordance with the 
half-integer winding number, and vanishes at = 0, indicating the presence of a single domain 
wall in the thin annulus. However, the condensate |'(/'(0)P = /^(0) presents a discontinuity 
in its angular derivative at the cut at = ivr, since we have 



f{ct>) = tanh^ ( , A /2(<^) = V2Kn,tanh ( \l - tanh^ ( ^ 



\V2 



(105) 



V\/2 / ' # \^/2 

Such a discontinuity is not physically acceptable, and is a consequence of the thin annulus 
limit which we have taken. Nevertheless, this approximate solution to the full LGH equations 
provides some interesting insight into the properties of half-integer vortex solutions. 

Let us now turn to the situation when the integration constant C does not vanish. In 
fact, the differential equation to be solved, (|101j ), is also that of a non relativistic particle of 
unit mass moving on the real line —00 < f < +00 and subjected to the potential 

Vif) = ^ + \^'ul ( f - \A , (106) 



2/2 ' 2 ''V 2 
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whose conserved energy Kj is thus 

Hence, there cannot exist a solution for f{(p) passing through the value / = for any finite 
value of the integration constant Kf, since the potential V{f) diverges at that point. An 
infinite value for Kf would rather be required, but this in turn is physically unacceptable 
since in such a case the angular variation df /d<j) of the order parameter would be infinite at 
all points except possibly at those where / = 0, thus also leading to an infinite energy. Hence 
we must conclude that no solution to the above equations may be found when C / 0. 

The physical interpretation of these results is as follows. Having taken the thin an- 
nulus limit Ua Ub with the restriction that j{(f)) = 0, this is consistent with the current 
conservation equation only for integer winding numbers without domain walls, since the axial 
symmetry of such solutions even for an annulus of finite width is such that j{u,<j)) = in 
any case. For half-integer winding numbers or integer vortices bound to domain walls how- 
ever, this limit is such that all states are squeezed out from the annulus, by acquiring an 
infinite energy. It is only when the applied field 6ext takes the distinguished values at which 
two consecutive integer winding states without domain walls become degenerate (which re- 
quires KUb > 1 /2) , that the states of half- integer winding number and with a single domain 
wall across the annulus retain a finite energy, but yet develop a discontinuity in the angular 
variation of the order parameter. 



4.3 Half-integer vortices in the disk 

Let us now finally turn to some exact solutions to the LGH equations (|6^ ) and (|65| ) in the disk, 
albeit numerical ones. The numerical resolution of these differential equations presents some 
genuine challenges. Even in the case of axially symmetric configurations, namely a single 
isolated giant vortex of integer winding number at the center of the disk, and thus in the 
absence of any domain wall, the direct numerical integration of the then ordinary differential 
equations for /(n), g{u) and b{u) — j{u) = in that case — are quite delicate. These equations 
then read 
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[l-f{n))f{u), (108) 



with of course ub{u) = dg{u)/du and the boundary conditions 
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, /(O) = ff L ^ 0, 



(109) 



L being the integer winding number. The main reason for difficulties in this case is, on the 
one hand, that the equations are non linear — the cubic term in in the LGH equation 
makes the stability of any numerical integration problematic — , and on the other hand, that 
the above local boundary conditions are in fact delocalized at the two boundaries of the disk. 
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namely at n = and at n = Uf,. Hence, whether integrating inwards or outwards, half of the 
initial values must be adjusted in order to satisfy the required boundary conditions at the 
other boundary. Nevertheless, it is still possible to manage such difficulties, and even obtain 
the annular vortex solutions although the precision requirements on the initial values then 
become critical|18]. 

In the case of configurations which are no longer axially symmetric, namely those in- 
volving either isolated integer vortices which are no longer at the center of the disk, or domain 
walls, or both, these difficulties are much greater, since the non linear equations are then par- 
tial differential ones in two variables for a larger number of quantities to be integrated, and 
with not only delocalized local boundary conditions but also the global boundary conditions 
corresponding to the integrated London equations for all contours in the plane. Furthermore, 
as we have seen, the current j{u, (p) is then also singular at the positions of vortices and inside 
the domain walls. These singularities at the location of a vortex which is at the center of the 
disk can be handled as above in the axially symmetric case, since we then have j{0,(j)) = 
and g{0, (p) = —Lq. Unfortunately, such a simplification does not apply to vortices away from 
the origin u = nor to domain walls. 

The alternative to integrating numerically the partial differential equations (^) and 
( |65[ ) subjected to the local and global boundary conditions, is to consider rather the numerical 
minimization of the free energy ( |73|) of the system using a steepest descent method of one 
type or another. Indeed, since the free energy coincides, up to a sign, with the action in 
the case of time independent configurations, the minimization of the energy is equivalent 
to solving the equations of motion. Most of the difficulties mentioned above may then be 
circumvented, by an appropriate choice of lattice discretization and of the positioning of the 
vortices and domain walls off the lattice sites. Such a procedure however, requires a lot of 
computational power, since all variables to be determined — in the present case the quantities 
j, g and / — are to be varied at each lattice site in order to follow the path of steepest descent 
in the free energy. This is, in a few words, the general method which has been applied in 



order to obtain the results to be presented here and in section 5^. Further details will be 
provided elsewhere. 

At this stage, simulations of only modest computer time have been run, leading to 
results which lack great precision. The only results to be described in this section apply for 
the following choice of parameters 

Ub = 3 , K = l. (110) 

The lattice discretization uses the polar coordinate parametrization with 10 intervals in the 
radial variable < u < Ub and 10 intervals in the angular range < (/> < tt. Indeed, the 
only configurations considered involve either a single integer vortex with L = 0, 1, 2 without 
a domain wall at the center of the disk, or else a single half- integer vortex with L = 1/2, 3/2 
also at the center of the disk and then bound onto only one domain wall. Furthermore, for 
reasons already given previously, we assume this domain wall to be straight and to extent 
all the way to the boundary of the disk by lying along one of its radii, which we take to be 
the one at (j) = 0. Consequently, all these configurations possess specific symmetry properties 
under (p —(p, which may be taken advantage of to reduce by half the calculation labor. 
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The rather coarse grained discretization used may imply a rather poor precision. Since 
the same procedure is apphed to all configurations, one may reasonably hope that whatever 
the numerical uncertainties, they would all contribute in the same direction so that at least 
the general properties, if not the precise numerical values, are to be trusted. Indeed, at least 
in the case of the configurations with integer winding number which are axially symmetric, the 
angular integration becomes trivial and a far more precise resolution, either by minimizing the 
free energy or by integrating the ordinary differential equations of motion using a 4th order 
Runge-Kutta method, becomes possible. The comparison between the different approaches 
then shows that the steepest descent method used in the half-disk leads to values for the 
free energy which are precise to 1% to 2% but which all deviate from their actual values in 
the same direction as expected. Hence, the same quality should apply to the results for the 
configurations with L = 1/2 and L = 3/2. 

Fig.l presents the results for the free energy of the disk in these different vortex and 
domain wall configurations as a function of the applied external magnetic field 6ext ■ Of course, 
the well known behaviour of the energy for the L = 0,1,2 configurations is reproduced. 
However, the curves for the solutions of winding numbers L = 1/2 and L = 3/2 are totally 
new, and are quite noteworthy. First, one may notice that at the values of 6ext where the 
curves of integer L cross each other, the energy values for the half-integer winding numbers 
always lie above those for the integer winding numbers, in agreement with the result of the 
analysis in the thin annulus limit. Second, the magnetic field dependency of the energy 
for half-integer winding numbers seems to be less pronounced than for the integer ones, 
and thus extends to larger values of 6ext before hitting the phase transition line ed, £ = 0. 
Nevertheless, these curves always lie above those of integer winding numbers (including L > 2 
not represented), so that the lowest energy states, at least for the geometry considered, always 
correspond to some giant vortex state without a domain wall. Third, the L = 1/2 and L = 3/2 
curves lie much closer to one another than do the integer winding number ones. In particular, 
it is quite remarkable that at 6ext = 0, both half-integer winding number curves lie below all 
the giant vortex ones without domain wall for L ^ [f^. 

Clearly, this is only one single example of such results in a very specific case, and it 
would be extremely interesting to see how these different features would change when varying 
the geometry of the disk and the value for k, when moving the vortices around within the 
disk, when allowing the domain walls to end on some other vortex inside the disk, when 
including more domain walls, and so on. Obviously, all these issues are of direct relevance 
to the detailed understanding of the magnetization properties and dynamics of mesoscopic 
superconducting disks and loops in varying electromagnetic fields. 

Even though not displayed, the numerical solutions for the configurations considered 
above also confirm the fact that there is indeed an electromagnetic current flow which quan- 
tum tunnels through the domain wall close to the vortex position for the L = 1/2,3/2 winding 
number cases. Furthermore, within the core of these vortices, the magnetic field is of course 
not axially symmetric, but in fact it even presents a shape such that its values first increase 
when moving out from the center of the vortex before decreasing towards the disk boundary. 

^''in this respect, it would be interesting to see where the curve for the L = 3/2 with three domain walls 
would lie. 
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This volcano-like shaped surface also possesses a valley-like small depression in the direction 
of the domain wall, along which the order parameter |/| takes of course smaller values and 
thus lets the magnetic field penetrate further in and decrease less rapidly as well (see (|65|)). 

In spite of the modesty of these first numerical results, they do establish the existence 
of half-integer winding number solutions to the LGH equations, bound to domain walls, and 
demonstrate the potential these new configurations have to contribute in crucial ways to the 
static and dynamic properties of any physical system which is described, even if only in an 
effective way, by the abelian U(l) Higgs model, and such instances of course are not confined 
to superconductors only. 



5 Bound 1/2-Domain Walls 
5.1 The solitonic limit 

Bound 1/2-domain walls were argued to provide the basic entities in terms of which any solu- 
tion to the LGH equations in the infinite plane could be viewed as a particular combination 
of such states being stacked together. Before considering in the next section the results of a 
modest attempt at a numerical resolution, here we shall present a specific bound 1/2-domain 



wall solution to the LGH equations ( |64D and (65) in the infinite plane. 



Even though this solution possesses an infinite energy, it may in fact be viewed to 
correspond to a bound 1/2-domain wall which is infinitely streched in a specific direction so 
that the two 1/2- vortices bound onto its edges are at infinity. It is this infinite length of the 
domain wall which explains its infinite energy of course. In spite of this fact however, the 
solution is interesting in that it illustrates some of the properties of bound 1/2-domain walls 
in a particularly simple limit, while also demonstrating that indeed the energy of domain 
walls is linear in their length when the vortices at their edges no longer overlap. 

To describe this specific solution, it is best to consider a system of cartesian coordinates 
X and y in the infinite plane (normalized to A), and to assume that the domain wall is aligned, 
say, with the x axis, so that the configuration is also invariant under translations parallel to 
that axis. Hence, all variables must be independent of x. Furthermore, given the physical 
interpretation, in which the 1/2-vortices are at infinity, one should expect that both the 
magnetic field b{y) as well as the current j{y), and thus also the electromagnetic current 
density J{y) = — f'^{y)j{y), vanish everywhere in the plane, 

h{y) = G , J(y) = 0. (Ill) 

These restrictions are clearly consistent with the current conservation equation, while the 
LGH equation for the order parameter then reduces to 

^/(y) = -K2(l-/2(y))/(y). (112) 

This latter equation is of course well known, and possesses solitonic solutions running from 
one vacuum / = — 1 of the Higgs potential to the other / = +1. In the present case, the 
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solution we are interested in is that corresponding to a single domain wall which separates 
two regions of the plane in which each of these two expectations values are reached at infinity. 

Hence, the infinitely streched bound f/2-domain wall is described by the one soliton or 
antisoliton configuration 

/(y) = ±tanh(^-^y). (113) 

Note how this solution takes values of opposite signs on both sides of the x axis, the choice 
of lb sign in this expression corresponding to a choice of sheet in the double sheeted covering 
of the plane which is defined by these one soliton and antisoliton configurations. Similarly, a 
> 1 soliton (or antisoliton) solution would correspond to N such infinitely streched bound 
1/2-domain walls all lying parallel to one another in the infinite plane. However, depending 
on the odd or even character of A^, such states are not topologically stable and would decay 
either to a single bound 1/2-domain wall if A'^ is odd, or else to the vacuum configuration 
|/(y)| = 1 if is even. 

Clearly, these configurations being translationally invariant along a specific direction in 
the plane, possess an energy ( [f3| ) which is infinite. Nevertheless, let us consider their energy 
per unit of length Ax along their symmetry axis, chosen here to be the x axis. One then 
finds for the one soliton configuration^ 




(114) 



Hence in this limit of an infinite separation of the two edges of the bound 1/2-domain wall, 
the energy does grow linearly with that distance. The amount of condensation energy — no 
magnetic energy density is present in the configuration — stored in the domain wall which is 
of constant thickness is proportional to its length. These are indeed the expected properties 
of streched bound 1/2-domain walls, as discussed in the Introduction. 



5.2 The bound 1/2-domain wall in the disk 

Let us now consider the numerical resolution of the LGH equations for a bound 1/2-domain 
wall in the disk topology, using the minimization of the free energy through steepest descent. 
Only the specific instance of a vanishing external magnetic field will be considered, 6ext = 0. 
Clearly, such an analysis touches directly onto the issue of the stability of the ANO vortex 
of winding number L = ±1, so that we shall in fact consider for different values of n the 
dependency of the energy (|73| ) of such states as a function of the separation between the two 
L = ±1/2 vortices bound onto the edges of such a domain wall. 

In the Introduction, it was argued that there might exist a critical value marking 
the boundary between the situation where, on the one hand, the bound 1/2-domain wall 
would be stable for a specific streched length (function of /t), and on the other hand, where 
it would always collapse into the L = ±1 ANO vortex. Furthermore, it was also suggested, 

^^The calculation is also feasable if the plane is of finite extent in the y direction, — j/o < y < yo, but the 
boundary conditions {df /dy){±yo) = are then not satisfied. 
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based on the interplay between repulsive and attractive magnetic and scalar forces and the 
situation in this respect for isolated giant vortices, that this critical value K\j2 could even 
coincide with the critical value k\ which governs a similar behaviour for integer vortices 
without domain walls, and which also coincides with the critical value Kc for which the latter 
states saturate the BPS bound. Based on such considerations, one would thus expect that 
the numerical resolution should possess some distinguishing feature when k = Kc = I/V2. In 
particular, it could then be the case that, for k < Kc the energy of the bound 1/2-domain wall 
would display a strictly positive monotonous increase as a function of its extension starting 
already from zero separation, while for k = Kc this same derivative would vanish, and finally 
for K > Kc, & small dip in the energy would be observed close to the origin corresponding to 
a split Abrikosov vortex at equilibrium, before turning over into the linear increase which is 
to be expected in any case for a sufficiently large separation. 

On the other hand, the existence of BPS lower bounds in the case of an infinite plane 
may also be brought to bear on this issue, in fact then leading to the opposite expectations! 
Indeed, even when k = Kc, we have seen that the only states which saturate the BPS bound 
are the isolated integer giant vortex ones, whereas any of the other vortices which involve at 
least one domain wall have an energy which is necessarily always strictly larger. Hence, in 
the infinite plane, the L = ±1 ANO vortex must be stable when k = Kc- For k > Kc however, 
since the BPS lower bound can never be saturated for whatever solution, the possibility 
that the energy could present a minimal value at non zero separation, before turning over 
into the linear increase, remains an open possibility. The same remark also applies when 
K < Kc, the more so since no BPS lower bound applies in this case, the term proportional to 
(k^ — 1/2) in ( [77| ) being then always strictly negative. In fact, from this point of view, this 
is the situation which would seem to offer the best case for an unstable ANO vortex since 
this term is the more negative the more the domain wall is streched (of course, the strictly 
positive two-dimensional volume integral in ([77| ) always ends up winning over for sufficient 
separation). These types of arguments would thus indicate that the value k = Kc could indeed 
play a distinguished role, but rather in exactly the opposite way to that which was argued 
for previously. 

To make things even more delicate, let us also note that for the more realistic case of 
a finite planar domain, such as a disk or an annulus, the BPS lower bounds no longer apply 
because of the contribution to the energy of the surface term at the boundary, so that even 
for the value the issue becomes uncertain. Furthermore, for finite domains, when 

the domain wall is sufficiently streched so that at least one of its edges comes close to the 
boundary, the Bean-Livingston barrier]^] would also become involved in the energy balance 
of the system, leading to further structure in the dependency of the energy on the vortex 
separation. 

Finally, let us also recall that in the singular limit k +00, all the solutions to the 
LGH equations become degenerate in energy and saturate the BPS bound. Since there are no 
interactions acting between the then membranes and Aharonov-Bohm fiux lines, the bound 
1/2-domain wall is stable for whatever separation of its two edges. This remark indicates 
that at least in this extreme regime, the ANO vortex may indeed split into such a state. 

Hence, one must conclude that given this lack of guidance for an educated guess, only a 
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detailed numerical analysis of the problem can answer this all important issue of the stability 
of the ANO vortex as a function of the coupling k. This section presents the results of a first 
attempt towards a resolution, which, unfortunately, proves to be still too modest to enable 
any definite conclusion. In fact, a full-fledged dedicated computationally intensive numerical 
analysis is required to resolve this important question, which we hope to be able to address 
elsewhere. 

The specific configuration which has been analyzed is as follows. Each of the L = 1/2 
vortices is placed at an equal distance d (in units of A) and on opposite sides from the center 
of the disk, while the straight domain wall joining them is aligned with the corresponding 
diameter, for the reasons given previously. Because of the four-fold symmetry of the configu- 
ration, it suffices to solve the problem in only one quadrant of the disk, say for < (f> < tt/2. 
The lattice discretization used 20 intervals in the radial direction < u < Ub and 16 intervals 
in the angular range < < 7r/2. Finally, to avoid as much as possible effects due to 
the finite boundary — such as the Bean-Livingston barrier — while still having a sufficiently 
fine-grained lattice discretization, the disk radius was set to 

Ub = 5. (115) 



The values for the coupling k which were considered are 

K : 0.5 , , 1.0 , 1.5. (116) 

For each of these values, the steepest descent minimization of the energy of the system was 
solved numerically, starting with d = 0, as a function of a regularly spaced scries of values 
for the distance d, which measures the position from the center of the disk of each of the two 
1/2- vortices bound onto the edges of the bound 1/2-domain wall thus streched to the length 
2d. 

Fig. 2 presents the results of these computations for the value 1/V2. The 

results for the other values of k show exactly the same behaviour — the only difference being 
in the vertical scale for energy values — , with in particular always a small dip in the energy 
at the first non zero step in the values for d whatever the value for k. Clearly, the appearance 
of such a dip for all the considered values of k does not agree with any of the expectations 
discussed above. It is not totally excluded that this may be due to subtle finite size effects 
because of the boundary at u = Ub, since = 6.74 • 10~^ is yet not that small. Indeed, the 
total magnetic flux ^[ub] was also monitored during the simulation and does also display a 
curious behaviour. When d = 0, one has a single L = I ANO vortex at the center of the disk, 
and ^[ub] then takes a value very close to unity as it should, namely ^[ub] = 0.989. However, 
when d then increases, the values for ^[ub] first increase getting even closer to the value of 
unity, before finally decreasing as they should when the two 1/2- vortices move closer to the 
disk boundary. On the other hand, it is not clear whether these two results for £ and ^[ub] 
could not also be a numerical artifact for small values of d, even though exactly the same 
computer code was applied to the situations with d = and d, ^ 0. Unfortunately at this 
stage, it is not possible to resolve that specific issue, which may only be addressed through 
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a full-fledged analysis including a much flner-grained discretization which is necessarily very 
much more computationally intensive. 

Nevertheless, note that the numerical solution docs display the expected linear increase 
in energy with separation, as soon as the bound 1/2-domain wall is sufficiently streched, 
namely when the two 1/2- vortices no longer overlap. Indeed, the linear behaviour sets in 
around d 2± 0.5 for all considered values of k (recall that d is measured in units of A). But 
here again, a much finer-grained analysis as a function of the separation is required. 

Hence, it seems fair at this stage to say that the fascinating and important issue of 
the possible instability of the ANO vortex as a function of k, which would then decay into 
a bound 1/2-domain wall streched to some length which should also depend on the value for 
K, is still completely open. Even though the quality of this first numerical analysis is too 
poor to reach a definite conclusion as well as a detailed understanding, the evidence which is 
available does not exclude the possibility, quite to the contrary. 

6 Conclusions and Outlook 

The fine-grained topological analysis which considers the transport of the complex scalar 
field around all possible finite contours in spacctimc, has thus uncovered a fascinatingly rich 
zoology of topological solutions to the Landau-Ginzburg-Higgs equations of the abelian U(l) 
Higgs model. Beyond the well known vortex states of integer winding number, the new 
configurations also include vortices of integer as well as half-integer winding number bound 
onto the edges of domain walls, all such vortex configurations possibly also being surrounded 
by annular current flows in the case of bounded spatial domains. The existence and physical 
consistency of these new states is related in a fundamental way to the U(l) local gauge 
invariance of the model, leading speciflcally to the all important fact that any of these states 
is to be characterized in terms of some double sheeted covering of the plane. The speciflc 
details of the physical properties of these solutions, and of their mathematical construction, 
have been described in the Introduction and established in the body of the paper. 

Even though the numerical analysis has not yet been able to answer some important 
issues, such as for instance the identification of the lowest energy configurations for specific 
combinations of such states, perhaps the most important of these cases being the possible 
instability of the usual Abrikosov-Nielsen-Olesen vortex to split into a domain wall with 
two vortices of winding number L = 1/2 bound onto its edges, the mere existence of such 
solutions raises a host of fascinating potential dynamical properties. Indeed, it now becomes 
possible to imagine that under specific electromagnetic disturbances, the magnetic vortices of 
systems described by these equations would not only move around in space, but could then 
also be deformed elastically into domain walls with vortices bound onto their edges, very 
much like rubber bands which would be pulled in different directions. Such properties should 
have important consequences, for instance, with regards to the transport of electromagnetic 
energy in the vacuum of systems obeying these equations, even if only in an effective way. 
Quite obviously, many issues of potential great direct physical interest, both in the static as 
well as in the dynamical regime, are raised by the existence of these new topological states 
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in the abelian U(l) Higgs model. 

Beyond their physical interest, these solutions raise other questions as well. From the 
mathematical physics point of view, it would presumably be very satisfying if the existence 
of such topological configurations, as well as of their stability, could be put on a firm mathe- 
matical ground, as has been done for the usual vortex solutions]^, One may also consider 
applying (an adapted form of) the fine-grained topological analysis to other field theories 
known to possess topological states, and possibly uncover further structure in these states 
as has been done here, for example by lifting specific symmetry properties such as rotational 
invariance (indeed, among the states described in this paper, only the usual vortices possess 
that symmetry). Extensions to higher spacetime dimensions or differential forms of higher 
degree could also be envisaged, having in mind the close analogy with D-branes and strings 
in M-theory, and the BPS bounds which are saturated for infinitely thin domain walls and 
vortices, corresponding then to sets of bound membranes and Aharonov-Bohm flux lines. 

Closer to particle physics considerations, these new solutions may also play an inte- 
resting role in the Landau-Ginzburg model of dual colour superconductivity as an effective 
description for the dual Meissner effect and confinement in QCD. Finally, since vortices and 
domain walls are quite generic in grand unified theories with spontaneous symmetry brea- 
king, the possibility that the states described here, or their cousins, could become involved 



in some interesting ways in texture formations! 25 1 in the Universe should not be excluded 
either. 

Nonetheless, the field for which the results of this paper offer perhaps the most im- 
mediate interest is that of superconductivity. The issue of the possible instability of the 
Abrikosov vortex is especially intriguing. Had we already been able to identify the equilibrium 
configuration of the streched bound 1/2-domain wall as a function of the Landau-Ginzburg 
parameter k, this could have been confronted to actual data, say in an Abrikosov lattice or in a 
nanoscopic superconducting diskf^. It may even be that the streching of the domain wall could 
be amplified to some degree by applying some external magnetic field. Such questions are 
certainly fascinating and deserve detailed study. Note that the recent results of Ref . [^] may 
be related to these issues. As remarked above, the existence of split vortices and domain walls, 
even if not as stable states, should also imply specific dynamical properties of superconductors 
submitted to electromagnetic disturbances, especially in regimes of frequencies whose wave 
lengths become comparable to the penetration and coherence lengths, and to the actual size 
of vortices and the thickness of domain walls on which electromagnetic waves would then pull 
almost as if rubber bands. 

In fact, we were led into the analysis of the issues of this paper because of our few years 
old project to develop a quantum detector for polarized particles based on superconducting 



loops 1 28 1 . Indeed, when such a loop is polarized in an external magnetic field close enough to 
the switching between two successive winding number states so that some local disturbance 
in the magnetic field would cause the loop to switch from one state to the other and back, 
this non linear response of the loop by a significant fraction of a full quantum of flux may 



^^A detailed magnetic field map of the bound 1/2-domain wall would also have to be known, in order to 
assess the required sensitivity in field values of the magnetic imaging technique to be used. Spatial resolution 
however, is no longer an issue |26[ 
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be picked up by an inductive circuit coupled to a SQUID. The same principle is at work 
in other devices using superconducting loops, albeit of larger dimensions, aimed towards 
the development of superconducting! 29 1 (or) quantum pO|] computers. However, in order to 
optimize the geometry of the loop with the particle detector application in mind, it becomes 
necessary to have a detailed understanding of the dynamics in a relativistic regime of the 
switching mechanism, which seems not be known. Only energy based considerations in the 
static regime are available for superconducting mesoscopic disks and annuli, mostly motivated 
by the results of Ref.[^. The only exception may be the analysis of Ref.|31| which provides 
a detailed picture for how an integer vortex may move in and out a superconducting disk 
or loop as a function of the applied field. Now that 1/2-vortices bound onto the edge of a 
1/2-domain wall have been uncovered from the same Landau-Ginzburg equations, the whole 
issue has to reconsidered again, to see whether there are specific loop geometries for which 
the sensitivity of the dynamical response to disturbances could be optimized, for example by 
having the L = 1/2 state have a lower energy than the L = 1 one when being degenerate 
with the L = state, so that magnetic flux would move in and out of such devices through 
a string of smaller bits than when only integer vortices are involved. 

Yet once again, the abelian U(l) Higgs model has proved to have still unknown precious 
tools hidden at the bottom of its treasure chest, to provide us, as if the Robinson Crusoes of 
the Universe stranded on our little "speckle of dust" , with the means to explore ever further 
the physical world at so many different scales. 
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Figure Captions 



Figure 1: The free energy E of the disk as a function of the apphed external magnetic field, 
for lift = 3 and k = 1. The vortex configurations are those discussed in section 4.3, with the 
vortices at the center of the disk and having winding numbers L = 0, 1, 2 and L = 1/2, 3/2, 
including in the latter two cases a single domain wall extending up to the disk boundary. 
The three dashed lines correspond, from bottom to top, to the states with L = 0, 1, 2 in the 
same order, while the two continuous lines correspond, in the same order again, to the states 
with L = 1/2, 3/2. For further details, see section f4.3| . 



Figure 2: The energy £ of the bound 1/2-domain wall in a disk with ui, = 5 and k = l/\/2 
as a function of the distance d to the disk center of each of the 1/2- vortices bound onto the 
edges of the domain wall. For further details, see section f).2\ The horizontal line is a guide 
to the eye, and corresponds to the value for d = 0. 
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